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We stndy the analytic properties of the scaling function associated with the 2D Ising model 
free energy in the critical domain T Tc, H —>■ 0. The analysis is based on numerical 
data obtained through the Truncated Free Fermion Space Approach. We determine the 
discontinuities across the Yang-Lee and Langer branch cuts. We conhrm the standard 
analyticity assumptions and propose “extended analyticity”; roughly speaking, the lat¬ 
ter states that the Yang-Lee branching point is the nearest singularity under hanger’s 
branch cut. We support the extended analyticity by evaluating numerically the associated 
“extended dispersion relation”. 
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1. Introduction 


The 2D Ising model is one of the best studied systems in statistical mechanics. 
Nonetheless, some questions concerning its criticality, notably in the presence of an exter¬ 
nal held H, remain open. The Ising model free energy exhibits a singularity at the critical 
point H = 0, T = Tc- The singularity is described in terms of the Euclidean quantum held 
theory known as the Ising Field Theory (IFT). It can be dehned as a perturbed conformal 
held theory through the action 

^ift = -^( 0 = 1 / 2 ) +'?■ J e{x)d^x + h j a{x)d‘^x, (1-1) 

where ^(c=i/ 2 ) stands for the action of c = 1/2 conformal held theory of free massless 
Majorana fermions, a{x) and e{x) are primary helds of conformal dimensions 1/16 and 1/2. 
To be precise, we assume that the normalizations of these helds are hxed by the usual CFT 
convention. 


(e(a:)e(0)) —> 1; ^ 1 


as 


\x\ 


0 . 


( 1 . 2 ) 


Under this normalization, the parameters r and h in (|1 . 1|) have mass dimensions 1 and 15/8, 
respectively. These parameters represent a deviation from the Ising model critical point. 


T = CrAT (1 + 0 (AT, if2)) ^ 
h = ChH (l + 0(AT,if2)) , 


(1.3a) 

(1.36) 


where AT = 1 — T/Tc, and real positive constants O,-, Ch (as well as the higher-order terms 
in the above relations) depend on the details of the microscopic (lattice) interactionS. The 
leading singular part Dsing{T, H) of the Ising model specihc free energy is universal, and 
it coincides with the vacuum energy density of the IFT (PD ; in what follows we use 
the notation T(27rr, h) for this quantity. It can be expressed through a universal scaling 
function $( 77 ) of a single variable - the scaling parameter 


ry = 27rr/h®/ 


15 


(1.4) 


(see Eq. (^.2CI| ) below for the precise dehnition of $( 77 )). This scaling function is of much 
interest as it controls all thermodynamic properties of the Ising model in the critical 


For the Ising models with nearest-neighbour interactions on square and triangular lattices, 
specific values of these constants can be found in for higher order terms in Eqs. (1.3) see 
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domain. Although there are many exact results (obtained through exact solutions of (Q) 
at = 0 and all r , and at r = 0 and all h P,|T3 ; these data are collected in |]^ ) 

as well as much numerical data about this function, its complete analytic 

characterization is still lacking. 

In this work we report preliminary results of a numerical study of the analytic 
properties of this scaling function. We use certain modification of the well-known Trun¬ 
cated Conformal Space Approach (TCSA) |T6,]T|, which we call Truncated Free-Fermion 
Space Approach (TFFSA). Although this modification is designed specificaly to treat the 
case of IFT (while TCSA is applicable to a wide class of perturbed CFT), in our case it 
produces better accuracy by taking full advantage of the fact that at /i = 0 the IFT ([Tt] ) 
is a free-fermion theory [^,|^. Using this approach, we compute numerically the scaling 
function $(? 7 ) for real as well as for some complex values of ry, to the accuracy sufficient to 
make exploration of its analytic structure in the complex 77 -plane. We locate the Yang-Lee 
edge singularity [^|^, and estimate some of its characteristics. Also, we study the free 
energy F'(27rr, h) at real r > 0 (low-temperature regime) and complex h] in particular, we 
determine, with reasonable accuracy, the imaginary part of the metastable branch of F 
for small as well as large values of h. For small h our result is in good agreement with 


the prediction from the critical droplet calculations moreover, we find a leading 

correction to the droplet model asymptotic. We formulate an “extended analyticity con¬ 
jecture” which states, roughly speaking, that the continuation of the free energy under 
hanger’s branch cut is analytic all the way down to the Yang-Lee singularity. We then 
use our numerical data to support this conjecture by verifying the associated “extended 
dispersion relation”. 

The paper is arranged as follows: In Sect. 2 some details of the TFFSA are presented. 
In Sect. 3 we discuss what is known and what is expected on the analytic properties of 
the free energy in the critical domain. Basic notations used throughout the paper are 
introduced there. The extended analyticity conjecture is formulated in Sect. 4, where the 
associated dispersion relation is also derived. In Sect. 5 we briefly discuss the excited 
(“meson”) states of IFT, and the role of the “false vacuum” resonance in the formation 
of the finite-size energy spectrum in this theory. The relation between the width of this 
resonance and the separation between the meson energy levels at the “near-intersection” 
points is presented there. Qualitative pictures of finite-size energy levels for real and pure 
imaginary h obtained through the TFFSA are described in Sect. 6, where we also explain 
in some detail how accurate numerical data for the scaling function are extracted from 
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the finite-size spectra. The data is described in Sect. 7, and in Sect. 8 we present our 
analytical interpretation of it. That includes approximations of the discontinuities across 
the Yang-Lee and Langer branch cuts, and numerical evaluation of the corresponding 
dispersion relations. Numerical support of the extended analyticity is also presented in 
this Section. In Sect. 9 we discuss possible physical significance of the extended analyticity. 


2. The Truncated Free-Fermion Space Approach 

As is well known (see e.g. 10). at zero external field the Ising model is equivalent to 
a free-fermion theory. Correspondingly, O can be written as 

•^IFT = *^FF + / cr^x) d^x , (2.1) 


Aff — 


1 


['tpd'ijj -)- 'ipd'ip -)- im d?x , m = 27rr . 


( 2 . 2 ) 


Here d = | — idy), d = \ -t- idy), where x — (x, y) are Cartesian coordinates, '0, ^ 

are chiral components of the Majorana fermi field, and a{x) is the “spin field” associated 
with this fermion. 


As in standard TCSA (see e.g. Ref. and references therein), we start with the 
IFT (0) in finite-size geometry, with one of the two Euclidean coordinates compactified 
on a circle of circumference i?, x-(-R ~ x. If y is treated as (Euclidean) time, the finite-size 
Hamiltonian associated with (p.l|) can be written as 




77ift = ^^FF + hV , V = 


a(x) dx , 


(2.3) 


where T^ff is the Hamiltonian of the free-fermion theory (|2.2| ). We are interested in the 
eigenvalues of T^ift, particularly in its ground-state energy Eq{R, m, fi), because for large R 
one expects to have 


Eq{R, m,h) = R E{m^ h) + O (exp (—Mi R)) , 


(2.4) 


with Ml being the gap in the spectrum of Mift aX R — oo, i.e. the mass of the lightest 


particle of the field theory (|1 . 1|) . In what follows, in refering to the eigenvalues of the 
Hamiltonian ( |2.3|) we typicaly use the notation E{R), with the arguments m, h suppressed; 


in particular, Eo{R) will stand for the ground-state eigenvalue of (|2.3|). 


3 







The free part H-pp of the Hamiltonian (^) is diagonal in the basis of iV-particle 
states of free fermions of mass \m\. At hnite R the space of states of (|2.2| ) splits into two 
sectors, the Neveu-Schwartz (NS) sector and Ramond (R) sector (with '0, ip antiperiodic 
or periodic as x —x + R, respectively). In each sector the particle momenta are quantized 
as Pn = 27: n/R where nG2 + l/2in NS sector, and n G 2 in R sector. In what follows 
we typicaly use the notation for integers, and ki for half-integers. The A^-particle states 
can be obtained from the NS and R vacua |0)ns and |0)r by applying the corresponding 
canonical fermionic creation operators. 


NS sector : |/ci, ■ • •,/cAr)NS = |0)ns /ci, • ■ ■,/cat G Z-h 1/2 , (2.5a) 

R sector : |ni, ■ • •, nAr)R = |0)r ni, ■ • •, G Z . (2.56) 


The normalizations of these states are hxed (up to phases) by conventional anticommuta¬ 
tors, 

{u/j, } ^k,k' 5 ^n,n' ■ (2-6) 

In all cases the energies associated with the A^-particle states, En^R), have standard form 

N 

(NS) (-R) = -^0 (NS) (-R) + ^ (R) , (2.7a) 

i=l 

N 

(R) (R) = Eq (R) E ^ ujrn (R) , (2-76) 

i=l 

where 

oJk(R) = \/Tn? (271 k/R)'^ ; oJn(R) = \/rrR + (27r n/R)^ (2-8) 

(with positive branch of the square root taken, in particular, c<;o(R) = \'kn\), and 

/ OO IQ . . 

— cosh6' log ^1 -t- , (2.9a) 

/ OO JCj . . 

— cosh6' log (1 - ^ ^2.96) 

are the eigenvalues associated with |0)ns and |0)r, respectively. The term 

9 

777/ 

R(m,0) = — logm^ (2-10) 

Stt 

in Eqs. (2.9) accounts for the famous Onsager’s singularity of the Ising free energy at 
zero h 0. 
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In order to treat the IFT with nonzero h we can admit only the states which respect the 
periodicity condition for the spin held, a{x + R,y) = +cr(x, y). This condition brings a dis¬ 
tinction between the cases m > 0 (the “high-T regime”) and m < 0 (the “low-T regime”). 
The admissible states are 


m > 0 : NS — states with N even , and R — states with N even , (2.11a) 

m < 0 : NS — states with N even , and R — states with N odd . (2.116) 


(at m = 0 the odd-A^ states in the R sector can be viewed as even-iV states, with one ao 
particle added). All the above statements are well-known (see e.g. [^ ). 

The operator h R in the full Hamiltonian ( p.3|) generates transitions between the states 
in NS and R sectors. Fortunately, all its matrix elements between the above states are 
known exactly. They are related in a simple way to the hnite-size formfactors of the 
held cr(x), for which an explicit expression exists. 


Ns(^i, ^ 2 , • • •, kKWiO, 0)|ni, n 2 , • • •, nAr)R 


K 


N 


S{R) n HOk,) n 9{0njFKM(^,^r--0kJ0. 


n\ 5 




( 2 . 12 ) 


i=i *=i 

where 6n {9k) stand for the hnite-size rapidities related to the integers n (half-integers k) 
by the equations 


|m|R sinh^fc = 27r/c, |m|R sinh^^ = 27rn . 


(2.13) 


In ( |2.12|) Fx,n is the well-known spin-held formfactor in inhnite-space |24 


Fk,n{9i, ■ ■ ■ ,9k\ 9[, - ■ ■, 9'j^) 


A 2 J (J 


tanh ^ ^ 

0<i<j<K 


where 


tanh 

0<p<q<N 




0<t<N 

(2.14) 


CT = s|m|'l®, s = = 1.35783834... (2.15) 

(with A standing for the Glaisher’s constant), and the rest of the factors represent hnite-size 
ehects. The overall factor S{R) is essentially the vacuum-vacuum matrix element 


aS{R) 


ns(0|(t(0, 0)|0)r form>0, 

Ns(0|ti(0,0)|0)R form<0. 


(2.16) 
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where ii{x) is the usual dual spin held |^, and S{R) is given by 
[mR)'^ fd9id92 sinh^i sinh 02 


S{R) = exp 


(27r)^ sinh(mi 2 cosh 0 i) sinh(mi?cosh 02 ) 


log 


coth 


9i - 92 


which was obtained in [^|. The momentum-dependent leg factors g and g are 


g{9) ~ /^/\m\R cosh 9 , g{9) = e /^^\m\R cosh 9 , 


where 


^ r 


1 


, —Irnli? cosh 9 ' 


2tt cosh( 6 * — 9') 


7T log 


1 + e 


— Imli? cosh 9 ' 


(2.17) 


(2.18) 


(2.19) 


The phase factors 2 ] where [...] denotes the integer part of the number) appearing 
in (|2.14|) can be removed by an appropriate phase rotation of the states, and thus play no 
role in the TFFSA computations. 

The above expression (p.l2| ) can be extracted from the result of recent papers . 
In fact, we have obtained it independently, before ||^ appeared, by a different approach. 
As our derivation seems to be simpler, we outline it in Appendix A. This expression is 


the m 7 ^ 0 generalization of corresponding massless matrix elements used in |T^ in the 
TCSA study of IFT with m = 0. 

It is useful to note that the Hamiltonian ( p.3|) can be rewritten to make its scaling 
form explicit, 

i^iFT = F^o(ns)(^) + \M Ho{r) + \m\^H^{r ), (2.20) 


where 




( 2 . 21 ) 


and the operators Ho{r) and Hfj{r) (corresponding to the terms H-p-p and hV in (|2.3|) ) 
depend on a dimensionless parameter r = \m\R only. 

The bulk energy density F{m,h) can be extracted from the large-i? asymptotic be¬ 
haviour (p.4|) of the ground-state eigenvalue of (p.3|) . As follows from ( p.2CI|) , it has the 
form 

, o o ^ / ..X , X 

( 2 . 22 ) 


F (m, h) = Ik— log -|- GR ), 
Stt 


where the hrst term is inherited from ( p.l0|) ; the scaling function G{^) (and related func¬ 
tion $( 77 ) dehned in ( p.20|) below) is the main object of our interest in this paper. 
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It is hardly possible to diagonalize the Hamiltonian (^) exactly. To render it tractable 
by nnmerical methods we follow the idea of TCSA, i.e. we nse hnite-dimensional approxi¬ 
mations (“trnncations”) of the inhnite dimensional space of the states (2.5), which inclnde 
only the states of snfhciently low energy. The trnncated Hamiltonian can be diagonalized 
nnmericaly, yielding an approximation to the gronnd-state energy f un ction Eo{R)-, the free 
energy is then extracted from its behavionr at snfhciently large R, according to ( |2.4|) . More 
details abont this analysis are presented in Sect. 6. 


3. Analytic properties of the scaling function 

In this section we discnss known results and expected analytic properties of the scaling 


function in (|2.22|) . 

In fact, to describe the free energy F(m, h) for both positive and negative m one needs 


two scaling functions G{^) in (|2.22|) . Although these two functions are analytically related 
(see Sect. 3.3 below) we use separate notations GiowiC) for m > 0 and GhighiO m < 0; 
these functions describe the free energy in low-T and high-T regimes, respectively. Both 
functions are dehned so that Giowi^) = Ghigh{0) = 0. 

3.1. The function GhighiO 

This function is even, i.e. GhighiO = Ghigh{—0- Around ^ = 0 it can be represented 
by a convergent power series in 

GhighiO = ^*2 + G 4 + Cg S - . (3.1) 



From 

References 

From high-T 
dispersion relation 

G 2 

-1.8452280... § 

-1.8452283 

G 4 

8.33370(1) 0 

8.33410 

Ge 

-95.1689(4) 0 

-95.1884 

Gs 

1457.55(11) 0 

1458.21 

Gio 

-25884(13) 0 

-25889 

Gi2 

5.03(1) X 10^ 0 

5.02 X 10^ 

Gi4 

— 

-1.04 X 10^ 


Table 1: Numerical values of the coefficients G 2 n in (^). The first column contains 
the data from Refs. [Q and |15]. The results obtained through high-T dispersion 
relation, Eq. (|3.5|) , with the use of our approximation (|8.2|) , are presented in the 
second column. 
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This series coincides with the perturbative expansion in h in the held theory (|1 . 1|) . 
The hrst coefficient G 2 is known through this perturbation theory exactly , and the coef- 
hcient G 4 was obtained with high accuracy by using formfactor expansion of the four-spin 
correlation function [^. There is a substantial amount of numerical data on the further 
coefficients. What appears to be rather accurate estimates of G 2 k up to G 12 are presented 
in a recent paper |^. The data on these coefficients are collected in Table 1. 

The function GhighiO can be analyticaly continued to complex values of The 
Yang-Lee theory [|T^[T^ guarantees analyticity of GhighiO the whole complex ^ plane 
with possible exception of the imaginary axis. At the imaginary axis one expects to 
observe branch cuts resulting from condensation of the Yang-Lee zeroes of the partition 
function in the thermodynamic limit. These branch cuts extend from i^o to ioo and from 
—i^o to —ioo, as shown in Fig. 1, where ^0 is some positive constant whose numerical value 
will be estimated in the Sect. 6, 

^0 = 0.18930(5). (3.2) 



Fig. 1: Analytic structure of GhighiO iu 
the complex ^ plane. YL and YL* denote 
the Yang-Lee edge singularities; associated 
branch cuts are shown as solid lines. 


Fig. 2 : Analytic structure of G;oii;(0- The 

solid line is the hanger’s branch cut. The 
star denotes the nearest singularity under 
this branch cut; its significance is discussed 
in Sect. 9. 


The branching points represent the Yang-Lee edge singularity |T^. Combining 
these analytic properties with the asymptotic behaviour GhighiO — (see Eq. (|3.34| ) 

below) one can derive the dispersion relation [IK 


GhighiO = 


2^mGimhit) ^ 

TT 


t{G+e) 


(3.3) 


$0 



















where the imaginary part of the function 


Gimhi,^') — it -|- 0 ) 


(3.4) 


relates in the usual way to the discontinuity across either of the branch cuts in Fig. 1. In 
particular, the coefficients G 2 n in (|3.1|) can be expressed in terms of this function as 


G2n = (-r 


Go 


2 Qm Gimh(t) ^ 

^2n+l jj- 


(3.5) 


In what follows we will refer to 


as the high-T dispersion relation. 


As is known ||^, the Yang-Lee edge singularity is a critical point, and the as¬ 
sociated CFT was identihed in |^] as the (nonunitary) minimal model with central 
charge cyl = —22/5. This CFT has only one relevant operator, the primary held of 
conformal dimension Ayl = —1/5. Therefore, one expects the singularities of GhighiO 
at ^ = ± z ^0 to be of the form 


Ghrskii) = Sa(A] + (fo + ^^) ” Sb{() + (?0 + “ SciO + subleading singular terms , 

(3.6) 

where the functions ^ regular at ^ = ±z^o- To some degree, 

the singular terms in ( |3.6|) can be understood in terms of the low-energy effective action 

Aflf =Aly CFT+ A(^) [ (j)_i/ 5 {x)(fx + [ {Tf){x)(fx 

J {2TimY J (3.7) 

-|- higher irrelevant operators , 


where (/>-i /5 is the relevant primary held mentioned above, and the rest of the action con¬ 
tains the contributions of irrelevant operators, the held (TT) (i.e. the L- 2 L -2 descendent 
of the identity operator) being the lowest of such operators. The dimensionless coupling 
constants A(^), q;(05 • • • 5 ^ire certain functions of the scaling parameter universal in the 
sense that they are uniquely determined by the original held theory (| 1 . 1 |) (but taken as 
“input” data in the ehective theory ( ^77|) ); all these functions are expected to be regular 
at ^ = ±z^o, and the critical point(s) ±z^o is dehned by the condition A(±z^o) = 0. Thus, 
in the vicinity of, say, — 

A(e = Ai («2 + ^7+o(({2 + ^")") ; 

*^(0 == <ao + o (^0 + ; 
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(3.8a) 

(3.86) 








and similarly for the other couplings. 

At ^ = —i^o (as well as at ^ = i^o) the held theory (|1.1| ) describes the Renormalization 
Group how from the Ising model hxed point with c = 1/2 down to the Yang-Lee hxed point 
with cyl = —22/5. In this case the mass gap Mi vanishes, and the hnite-size energy levels 
are expected to approach the linear asymptotic (|2.4| ) with a power-like, not exponential, 
accuracy. These power-like corrections to the linear asymptotic (|2.4|) can be analysed 
through the perturbation theory based on the ehective action (|3.7| ). In particular, for the 
ground-state level of (|2.20|) with ^ one obtains 


o{R) ^®^24(mR)2 (24)2 (mR)4 




(3.9) 


where Cefi = cyl — 24 Ayl = 2/5 is the “ehective central charge” of the Yang-Lee CFT, 
CKO is the leading coefficient in (|3.8d), and Fq stands for the free energy Fq associated with 


this how; the last quantity is related to the value of the scaling function GhighiO 
singular point, 

Fq = F{m,0)+m'^ Ghighi±i^o) ■ (3.10) 


The calculation leading to 


is similar to that presented in [31|; we skip it in this report. 


3.2. The function GiowiO 

Unlike GhighiO above, the analytic continuation of GiowiO ^ot an even function 
of More precisely, for m > 0 the analytic continuation in h of the free energy (|2.22| ) 
yields two diherent analytic functions, depending on whether one continues from positive or 
negative parts of the real axis. In what follows, by GiowiO always understand analytic 
continuation of the free energy from positive part of the real ^ axis. Again, according to the 
Yang-Lee theory, GiowiC) is analytic in the right half-plane 3fJe^ >0. hanger’s theory m 
(as well as earlier calculations within “droplet models” |^,0) predicts a weak singularity 


at ^ = 0, and it is usually assumed (see e.g. [p^ , |35| , p6| ) that GiowiO i^ analytic in the full 
complex ^-plane with the branch cut from 0 to —oo, as is shown in Fig. 2. We will call this 
the standard analyticity assumption. This assumption will be conhrmed by our numerical 
analysis in Sect. 7. The function GiowiO admits an asymptotic expansion in powers of 


Glow (0 — Gi ^ + G 2 + 


(3.11) 
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Again, the coefficients Gn here in principle can be computed by means of the perturbation 


theory in h in (|2.1|) . Thus, the first coefficient is directly related to the spontaneous 
magnetization at zero h, given by (|2.15|), 


Gi = -s. 


(3.12) 


The coefficient G 2 is also known exactly, since its computation involves integrating the 
two-spin correlation function, which is determined in terms of the Painleve functions R]. 


Several further coefficients were estimated in using exact large-distance expansions of 
multi-spin correlation functions at zero h. Numerical estimates from lattice series analysis 
are also available 0 . We quote these results in Table 2. 



From 

References 

From low-T 

dispersion relation 

Gi 

-1.35783834... 

-1.35783835 

G 2 

-0.0489532... 0 

-0.0489589 

Gs 

0.03875291^; 0.039(1) ^ 

0.0388954 

G 4 

-0.0685535 0; -0.0685(2)® 

-0.0685060 

G 5 

— 

0.18453 

Ge 

— 

-0.66215 

G 7 

— 

2.952 

Gs 

— 

-15.69 

Gg 

— 

96.76 

Gio 

— 

-6.79 X 102 

Gn 

— 

5.34 X lO^ 

Gi2 

— 

-4.66 X 10^ 

Gi3 

— 

4.46 X 10® 

Gi4 

— 

-4.66 X 10® 


Table 2: Numerical values of the coefficients in (3.11). Exact results for 
Gi and G 2 , and nnmerical estimates fonnd in literatnre, are collected in the first 
column. The second column shows resnlts obtained throngh low-T dispersion rela¬ 
tion, Eq. ( 3.13 ), using onr approximation ( 3.14 ). 


If one makes the above standard analyticity assumption (as in Fig. 2), the coeffi¬ 
cients Gn can be represented as the integrals 

n+l Am Gmetajt) ^ 

./n TT 


G. = (-) 


n = 2, 3,4,..., 


(3.13) 
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where the function 


Gmeta{(,) = Glow{—^ + * 0 ) 


(3.14) 


at positive ^ describes the values of GiowiO ^^al negative at the upper edge of the 
branch cut in Fig. 2. Correspondingly, the discontinuity of GiowiO across this branch cut 
equals 2i 9m Gmeta{—0- The Eq. (|3.13|) follows directly from the dispersion relation 

^mGmetait) dt 


GioU0=Gi^-^^ 


{t + 0 ^ ’ 


(3.15) 


where the integral converges at any hnite ^ because GiowiO ~ as ^ ^ oo (as follows 

from (|3.28|) and ( |3.29|) below). 

The function GmetaiO deservingly attracts much attention, for at least two reasons. 
One is the commonly accepted interpretation of the function 


Fmetaij^^ F) — TTl | log m T G^netai^) 




(3.16) 


which coincides with the analytic continuation of the free energy (|2.22|) to negative values 
of /i, as the free energy associated with the metastable state at T < T^. In fact, to the 
best of our knowledge, it is the only mathematically precise dehnition of the metastable 
free energy available today. The function ( 3.16| ) takes complex values, and its imaginary 
part 

T{m,h) = Fmetaim, h) = rri^ ^mGmetaiO (3-17) 

is interpreted as the rate of decay of the metastable state. More precisely, according to 
Langer’s theory the actual rate of decay is asymptotically proportional (with the 

constant factor absorbing the time scale) to this imaginary part in the limit h —> 0 (or 
^ > 0). It is tempting to assume (as is often done in the literature on nucleation theory) 

that similar relation extends to some hnite domain of h, although it is understood that 
the proportionality coefficient, being sensitive to at least some details of the kinetic model, 
may very well depend on h, and its degree of universality is not clear. Another, and much 
better understood, interpretation of ( 3.16|) is in terms of the quantum held theory (|1.1J ) 
in 1 + 1 Minkowski space-time. Namely, at r > 0 and h ^ 0 this held theory exhibits 
a global resonance state commonly refered to as the “false vacuum”, and the quantity 
( P.16|) coincides with the associated (complex) energy density; in particular, the imaginary 


part of ( |3.16|) gives precisely the decay probability (per unit volume and per unit time) 
of this resonance state |]39|,|^,^ . If h is small, both of the above interpretations allow 
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one to justify the validity of the instanton saddle-point calculation, which in the D = 2 
case yields the following h —>■ 0 asymptotic behaviour of the imaginary part of ( 3.16 ) [pT 
(see |p^ , ^ ^ ,|4T| , ^ for such calculations in more general context), 


r(m, h) 


a h 


exp 


TT m 

2 ah 


as /i ^ 0 , 


(3.18) 


where a = s is the magnetization at zero h (see Eq. ( p.l5|) ). It is fair to say that there 

is still some controversy about whether ( p.l 8 |) gives the correct numerical coefficient in the 
asymptotic behaviour of this imaginary part [^. Our analysis in Sect. 7 is completely 
consistent with the coefficient in ( p.l 8 |) , and actually offers the leading correction to this 
asymptotic (see Eq. (^T^)). 


3.3. The function $( 77 ) 

Although the scaling functions Ghigh and Glow above are dehned independently, in 
fact they can be analytically related one to another. To do this it is useful to introduce 
another scaling variable, 

?7 = , (3.19) 

and rewrite the free energy (| 2 . 22 D as 


m 


F(m, h) = — logm^ -|- 
Stt 


16/15 $(^)_ I'g ^O) 

Of course, the scaling function $( 77 ) here is related to the function(s) G{ff) in (|2.22|) , and 
vice versa. If both m and h are real and positive we have 

V =, 

and therefore 

for real 77 > 0 . 


$( 77 ) = r]‘^Giow 

Similarly, if h > 0 but m < 0 the variables (| 2 . 21 | ) and (|3.19|) are related as 

so that 

for real 77 < 0 . 


(3.21) 

(3.22) 


(3.23) 

(3.24) 


$(?7) = ri'^ Ghigh (l/(-77)^^/®) 

On the other hand, for hxed h ^ 0 the free energy F{m, h), as function of m, is expected 

to be analytic at all hnite real tti, including the point tu = 0 , for if h 7 ^ 0 the correlation 

length remains hnite even at tu = 0. Therefore the scaling function $( 77 ) can be written 

as 2 

^(17) =+ ^(17): ( 3 - 25 ) 

Stt 

where the fu n ction $( 77 ) is analytic at all hnite real 77 . 
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Fig. 3: Mapping of Fig. 1 under the vari¬ 
able transformation ( 3.23D ; the wedge HTW 
(—< arg(—ry) < is the image of the 
right half-plane in Fig. 1. Rotating the branch 
cuts as suggested by the arrows opens the prin¬ 
cipal sheet of the ry-plane shown in Fig. 4. 


Fig. 4: Principal sheet of the 77 -plane for 
the scaling function ^( 77 ). The wedge LTW 
(—< arg( 77 ) < is the image of the 
principal sheet in Fig. 2 under the map ( 3.21 ). 
Analyticity of ^( 77 ) in the wedge SHD 
< arg( 77 ) < — is the subject of the 
’’extended analyticity” conjecture in Sect. 4. 


Eqs. (|3.22|) and ( |3.24|) can be promoted to analytic relations valid for complex valnes 
of 77 . One notes that the variable transformation (|3.23|) maps the right half-plane in Fig. 1 
onto the wedge —47r/15 < arg(— 77 ) < 47r/15; we call it the “High-T Wedge”, see Fig. 3, 
where it is marked as the region HTW. Therefore the following analytic relation holds in 
this wedge, 


2 

-^\og{-rii) + ^{ri) = ifGhigh(^/{-vY^^^^ for - 47r/15 < arg(-? 7 ) < 47r/15 , 

(3.26) 

where the principle branch of the log fu n ction (i.e. log( 77 ) is real for real 77 > 0 ) shonld be 
taken. The Yang-Lee branching points are located at the points 77 = — in the 

complex 77 -plane, where 

Yo = ■ (3-27) 


The way the corresponding branch cuts are drawn in Fig. 1 corresponds to placing 
them along the rays 77 = with y running from —Yq to -|-cx), as is done in Fig. 3. 

However, it is convenient to rotate these branch cuts 180°, as is suggested by arrows 
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in Fig. 3. This way one exposes the fnll “principal sheet” of the complex ry plane, shown 
in Fig. 4, where the branch cuts extend along the rays ry = —oo < y < —Yq. 

In this picture the way to reach the left half plane in Fig. 1 is to go from the region 
HTW under either of the cuts. This complication is not important, as in view of the 
symmetry GhighiO = Ghigh{—0 expects to see there merely another copy of the 

principal sheet. More importantly, there is now a direct way (through some vicinity of the 
real axis, where 4>(?y) is analytic) into the right-hand part of the principal sheet. Here, in 
the wedge —47r/15 < arg(?y) < 47r/15 (the domain LTW_|_), which is the image of the right 
half plane in Fig. 2 under the variable transformation ( p.21|) , the function 4>(7y) is related 
to GiowiO as 
2 

— ^ log(7y) -I- l'(7y) = 7y^ Gzoiy for — 47r/15 < arg(7y) < 47r/15 , (3.28) 

where again the principle branch of the logarithm is understood. Thus, 4>(?y) is analytic 
in the wedges HTW and LTW+ in Fig. 4, where it matches the functions Ghigh and Glow 
through (|3.26[ ) and (|3.28| ), plus it is analytic in some hnite neighbourhood of the 
point 1 ] = 0. Moreover, the transformation (|3.21|) maps the full principal sheet of the 
complex ^-plane in Fig. 2 onto twice wider “Low-T Wedge” —Stt/IS < arg(?y) < Stt/IS 
(shown as LTW in Fig. 4), with the upper and lower edges of the branch cut in Fig. 2 
represented by the rays arg(?y) = =f Stt/IS. The standard analyticity assumption (see 
Sect. 3.2) implies analyticity of 4’(Ty) in this wider wedge. Not unexpectedly, our numerical 
analysis in Sect. 8 will conhrm this to an high degree of precision. The question interesting 
indeed is about analytic properties of this scaling f un ction in the remaining part of the 
7y-plane, the “Shadow Domain” —IItt/IS < arg(?y) < —Stt/IS (SHD in Fig. 4). Signihcance 
of this problem will be discussed in Sect. 9. In Sect. 4 below we will formulate what we 
believe to be the most natural conjecture about the shadow domain analyticity. 

In any case, 4’(ry) is analytic at ?y = 0, and in some hnite domain around this point it 
can be represented as a convergent power series, 

4'(ry) = $0 + $iry$ 2 ry^ 4-. (3.29) 

In fact, the hrst two coefficients in series (|3.29|) are known exactly, thanks to integrability 
of the held theory (|1 . 1|) with r = 0 |p . The coefficient 4>o appears as the amplitude in its 
free energy, 

F(0,h) = $0 (3.30) 


IS 











which was computed in [10 


r(l/3)r(l/5)r(7/15) / 4xU^(13/16)r(3/4) y^‘» 

° 27rr(2/3)r(4/5)r(8/15) V r2(3/16)r(l/4) J ■ ^ ^ 

The coefficient $i is related to the expectation value of the energy density field e{x) in 
this integrable theory, 

(e)^=o = -27r$l|/l|"/l^ (3.32) 


which was determined in [^ 


(e) 1^=0 = (2.00314...) . (3.33) 

Further coefficients can in principle be computed using perturbation theory in r; this 
however is a hard problem because multipoint correlation functions of the r = 0 theory 
are not known exactly (see however Ref. [^). Our numerical analysis in Sect. 6 yields 
next six coefficients $2 to $7 with reasonable accuracy (see Table 3). 



From 

References 

From 

TFFSA data 

From extended 
dispersion relation 

$0 

-1.1977334... [N 

-1.1977331 

-1.1977320 

$1 

-0.3188096... til 

-0.3188103 

-0.3188192 

$2 

— 

0.1108867 

0.1108915 

$3 

— 

0.0164266 

0.0164252 

$4 

— 

-2.64 X 10“^ 

-2.64 X 10“'^ 

$5 

— 

-5.14 X 10“^ 

-5.14 X 10-4 

$6 

— 

2.07 X 10“^ 

2.09 X 10-4 

$7 

— 

-4.52 X 10“^ 

-4.48 X 10-^ 

$8 

— 

— 

3.16 X 10-^ 

$9 

— 

— 

4.31 X 10-6 

$10 

— 

— 

-1.99 X 10-6 


Table 3: First few coefficients of the power series (3.29): Exact results from in¬ 
tegrable field theory (first column), estimates from TFFSA data (second column), 
and results obtained from the extended dispersion relation, Eq. (4.10), using ap¬ 
proximation (8.21|) (third column). 
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According to (|3.26|) and (|3.28|) , the coefficients control the asymptotics of the 
functions GhighiO GiowiO large for instance, for —7r/2 < arg(^) < 7 r /2 and 
sufficiently large |^| 


GhighiO 


8C2-n) 

15 as ^ > cxo . (3.34) 


n=2 


Similar expansion (without the factors) holds for GiowiO ^^e wider domain 

—71 < arg(^) < 71. Conversely, expansions (|3.1|) and (|3.11|) can be rewritten as large -?7 ex¬ 
pansions of the scaling f un ction $( 77 ) in the corresponding domains of the 77 plane. Thus, 
the expansion 


47r 


log(-77) + $( 77 ) = G 2 {-rj) 4 Gi (- 77 ) 4 Gg {-vY 


37 


„ 47r 47r 

tor-< argi— 77 ) < — 

15 ^ 15 


(3.35) 


converges to $( 77 ) for sufficiently large 77 in the domain HTW in Fig. 4. In the domain 
LTW the function $( 77 ) has asymptotic expansion 


77 ~ ~1~7~29 

-log( 77 ) -|- $(77) ~ Gi 778 -|- G2 77“4 -|- (^3 -|- ... as rj ^ 00 , 

Att 


for 


Stt 

15 


Stt 

<arg(77)<—. 

15 


(3.36) 


To prepare notations for our analysis in Sects. 6-7 let us briefly describe the expected 
behaviour of the scaling function $( 77 ) along the axes AOB*, A*OB and on the rays OC, 
OC* in Fig. 4; these domains are of particular interest for our discussion. 

The axis AOB* corresponds to 77 = 6^*77 with real (positive as well as negative) y. 
The values of the scaling function $( 77 ) along this axis are related to the free energy F{m, h) 
at pure imaginary h as follows. Let us introduce the function 


^^mh{y) = -|^ log \y\ - ^ 


+ A 




(3.37) 


According to ( |2.22|) , (|3.2CI|) and ( p.25|) , this function relates to the analytic continuations 
of GhighiO and GiowiO to pure imaginary Namely, 


^imh{y) = y‘^ Glow + oj for 77 > 0 ; 

^imhiy) = y^ Ghigh + o) 
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for 77 < 0 . 


(3.38a) 

(3.386) 


















Therefore it takes real values for y < —Yq (the upper edge of the branch cut A-YL in 
Fig. 4), and becomes complex-valued for y > —Yq, where its imaginary part relates to 
imaginary part of the function (|3.4|), 


(3.39) 


Dimhiy) = ^rn^imhiy) = y'^^rnGimh y<0. 


The structure of the singularity of ^imh{y) at —Yq can be described in terms of the 


expansion ( p.6|) . As in the analysis in Sects. 6-8 we use the variable ry rather then let 
us introduce here more suitable notations. In the vicinity of —Yq we expect to have 


^imh{y) = A{y) + B{y) (-Yq - y - tO) e C{y) (-Yq - y - tO) 3 + 


(3.40) 


where the real analytic functions A(y), B{y), C{y), ..., are regular at y = — Yq, i.e. they 
admit the power-series expansions 


A{y) = Ao + Ai (-Yo - y) + A 2 (-Yq - yf -\ -, 

B{y) = Bq + Bi (—Yq — y) + ■ ■ ■ , 

C{y)^Co + ---, 


(3.41a) 

(3.416) 

(3.41c) 


around the point — Yq, with real coefficients. The structure ( p.40|) follows from the expected 
form (|3.7| ) of the effective action in the vicinity of the critical point. The hrst subleading 
singular term written down in (|3.40|) is due to the term TT in (|3.7D . Simple perturbative 
analysis (which we skip in this report) allows one to relate some of the coefficients in (3.41) 
to the parameters in (3.8), 


Bq — /yl o^yl (15Ai Y(^ 


n5/6 


C'o = aoSoV4n^ 


(3.42) 


where /yl = —"s/S/lS is the free energy amplitude of the integrable Yang-Lee held 
theory P5|,|46[|, and uyl = 2.6429446... is the corresponding amplitude in its “mass to 


coupling relation” [46.47 


Let ns also write down the y-series expansions of (|3.39|) , which follow directly 


from (|3.29|) . For sufficiently small |y| we have 


Dimh{y) — ^ ^ ^ny 

n=0 


(3.43) 
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where 


for n 7 ^ 2 ; 


(3.44) 


5 ^ . 47r (n - 2) 

sin 


$2 = 


1 


15 ' ’ 15 

The axis OC in Fig. 4 (i.e. tj = y with real y > 0) is the image of the upper edge 
of the branch cut in Fig. 2. The values of $(? 7 ) along this axis are related to the complex 
free energy (|3.16| ), 

2iy‘^ 


imetaiy) = “ |^ '“S S' + I* ( S'" ‘'* ) = 'l/G-^eta ( l/v' 


(3.45) 


In view of (|3.15| ) and (|3.17|) , the imaginary part of this fu n ction plays particularly important 
role, so we introduce for it a separate notation 


Dmetaiy) = ^metaiv) = ^ ^ ^ (V ^ j - 6 * $ ( y e* 

As follows from (|3.45|) , this f un ction enjoys a power series expansion 

CXD 

Dmetaiy) = Dny'^ , 

with 


(3.46) 


n=0 


Dr,. = sin 


Stt (2 — n) 


for n 7 ^ 2 ; D 2 = 


(3.47) 


(3.48) 


15 ' 15 

which converges for sufficiently small y. On the other hand, Dmetaiy) decays very fast at 
large positive y. According to ( p.l8|) . 


Dmetaiy) ^ +00 , 

271 


(3.49) 


where Viy) stands for the function 

Viy) = y^ exp (^-7 ry^/2sj 


(3.50) 


4. Extended Analyticity 

In the previous Section we have described the domains of analyticity of the scaling 
function $( 77 ) which follow from previously known results. An unexplored area lays between 
the High-T and Low-T Wedges in Fig. 4, the “Shadow Domain” SHD. In more conventional 
notations (i.e. in terms of the scaling variable ^) this domain is located under the banger’s 
branch cut in Fig. 2, imediately below the negative real axis. Alternatively, it can be 
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reached by going under either of the branch cuts in Fig. 1. We believe the problem of 
establishing the analytic properties of free energy in the shadow domain is an important 
one, and it was one of the motivations of this work. Some reasons for our interest in this 
problem are explained in Sect. 9. 

The simplest conceivable possibility is that the f un ction $( 77 ) is actually analytic in 
the wedge —IItt/IS < arg(? 7 ) < ll7r/15 which includes the whole of the shadow domain. 
It is also the most elegant possibility (to our sense of beauty, that is), because it renders 
the Yang-Lee point an additional signihcance of being the closest singularity under the 
banger’s branch cut in Fig. 2. We formulate it as the following 

Extended Analyticity Conjecture: The scaling function $( 77 ) is analytic in the whole 
complex ? 7 -plane with two branch cuts extending from the points YL and YL* to inhnity, 
along the rays YL — A and YL* — A* in Fig. 4, respectively. 

The extended analyticity is equivalent to certain “extended” dispersion relation, ex¬ 
pressing this scaling function in terms of its discontinuity across the branch cuts in Fig. 4. 
To write it down, let us introduce the f un ction 

A{y) = i (- 7 / - $ (^ 7 / ^ , (4.1) 

which is the discontinuity of 4 ’(r 7 ) across the branch cut YL — A, with additional phase 
factor introduced for later convenience. Here 77 is a positive real variable interpreted as the 
coordinate along the ray O — A in Fig. 4. The fu n ction ( [4.1|) vanishes for y < Yq, while 
at y > Yq it takes complex values. In the last domain it coincides with certain analytic 
continuation of the function Dimh{y), 

A( Yo + ^) = D,mh {-Yo + ^ , (z > 0 ). (4.2) 

Unfortunately, this function grows too fast at large positive y {^{y) —y^/4 as 

y —>■ +cxd), and writing down the dispersion integral directly for ^{rj) would require at 
least three subtractions. It is more convenient to use instead the f un ction 

A{y) = Aiog ( 77 ) + A( 77 ) , (4.3) 

where 

A;og(77) = ^77^. (4.4) 
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The function (O) is related to the mismatch between the functions 


and 


ihish(tl) = -^ log(-i)) + 


ilow(ri) =-^ tog(rj) + $( 17 ) 


(4.5) 


(4.6) 


at the ray O — A. According to the asymptotics (|3.35|) and (|3.36|) , the function 
exhibits a much slower growth rate, A(y) Aass(y) as ?/ —> +cx), where 


Aass{y) = -e *8 Giys , 

and hence one can write down the dispersion integral with only one subtraction, 

^ roo y^e(^e^ A{y)^+r]^e(^e^ A{y)^ 


(4.7) 


$( 77 ) = $0 + ^logiv) - 


TT 


'Vo 


7/2 + 2yr/ cos(47r/15) + 


where the term 

^log(y) = 




cos(47r/15) + — log (Yq + 2 yYo cos(47r/15) + rf) 

47 r OTT 


(4.8) 


(4.9) 


comes from the domain 0 < ?/ < To ia which A{y) coincides with (^4.4|) . 

We will refer to Eq. (^4.8| ) as the “extended dispersion relation”. Its validity is 
equivalent to the validity of the extended analyticity assumption stated above, and under 
this assumption Eq. ( |I.8|) must hold in the whole complex ry-plane, including the shadow 
domain SHD in Fig. 4. 

From ( [4.8|) one readily derives the following expressions for the coefficients of the 
Taylor expansion (|3.29|) , 

A{y)-Aass{y) 


$0 = 3fJe 


Stt^ 

e 15 


TT 


y 


dy 


$1 = 


.Xy 

4ti 


cos(47r/15) — 3fJe 


f- 47 ri 

e 15 


A{y) dy 


TT 


'Yo 


$2 = 1 — log To + 
47r 


7*00 


A{y)dy 




^(2-n) 


e 15 


'Yo y 

r (A(y) 


Aiog{y))dy 


TT 


,n+l 


for n > 2 . 


(4.10a) 

(4.106) 

(4.10c) 

(4.10d) 


Jyo y 

In Sect. 8 below we will bring up some numerical evidence to the validity of this 
extended dispersion relation, thus furnishing a support for the above extended analyticity 
conjecture. 
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5. Excited states. False vacuum 


Although in this work we are mostly interested in the free energy F{m, h), and hence 
concentrate most of attention on the ground-state energy (|2.4|) of the finite-size Hamil¬ 


tonian (|2.3|) , some analysis of excited states will prove to be useful in Sect. 8, where we 
study the function ^metaiy)- To prepare for this discussion let us briefly remind qualita¬ 
tive properties of the quantum field theory (|2.1J ), in infinite space-time, interpreted as a 


particle theory [48 


For m > 0 and h = 0 (i.e. for ry = -|-cxo) we have just a free fermion theory 
In this case the vacuum of the bulk theory is two-fold degenerate (the ground states in 
NS and R sectors degenerate as R —> cxo), thus manifesting the spontaneous breakdown 
of the a —> —a symmetry. The particles (i.e. the free fermions of (|2.2|) ) are identified 
with the “domain walls”. Adding the interaction term in (|2.1J ) creates a confining force 
between these fermions (which henceforth are refered to as “quarks”), and for nonzero h 
the particle spectrum contains only their bound states - “mesons”. If m > 0 and h is small 
(i.e. ry is large positive) the corresponding “string tension” is small (~ 2 ah, where a is 
the spontaneous magnetization at zero field given by (p.l5| )), and there is a large number 
of stable mesons, their masses Mi densely filling the interval between 2m and 4m (the 
mesons with M > 2 Mi, Mi being the lightest meson mass, are generally unstable). In 
this region of small h the lower part of the meson mass spectrum (with Mi — 2m m) 
can be understood in terms of non-relativistic quarks interacting via a linear conhning 


potential. This interpretation agrees with the result of for the masses of these lightest 
mesons, 

2 

(2 a h)3 Zi , , , 

Mi —2 m—>•-j- as —> 0 , (5.1) 


m3 


where —Zi, i — 1,2,3,..., are zeroes of the Airy function, Ki(—Zi) = 0. In fact, a few hrst 
relativistic corrections to these masses can be computed (see Appendix B), 


Mi = m <2 -\- 


(2s) 3 Zi 

5 

r/4 


(2s)^zi fllz! 

V 1400 


20?y2 




(5.2) 


where the constant q 2 is the same as in Eq. ( b.8|) below. 

As rj decreases, the heavier of the mesons gradualy disappear from the spectrum of 
stable particles; when their masses exceed the stability threshold 2 Mi they (in general) 
become resonance states. In particular, when rj reaches zero (i.e. m = 0), only three 
stable mesons remain below the threshold; in this case however there are additional hve 
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stable particles above the threshold, which owe their stability to the fact that the held 
theory (O) with r = 0 is integrable 0 . This process continues when ry becomes negative, 
until hnally at ry < 772 (?y 2 ~ —2.09, see Sect. 6 below) only one particle remains stable. As 
ij —>■ —00 its mass Mi approaches |m| i. 

Ml = |m| ^1 + a/(—?y)^ + O ^(—ry)“^ j ^ ; a 10.75 . (5.3) 

Although a detailed discussion of the mass spectrum is outside the scope of this paper (we 
intend to present it separately), we show in Fig. 5 the rj dependence of the hrst few meson 
masses obtained using the TFFSA. 
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Fig. 5 : Masses of the three lighest particles in the Ising field theory (O). The solid 
lines are the plots of the dimensionless ratios Mi/\h\^^^^ (i = 1, 2, 3) versus the pa¬ 
rameter ( b-isp . The dashed lines represent the corresponding large-|ry| expansions, 
Eqs. (b-2D and (b-3D (with all terms explicitely written in these Eqs. included). 


^ The numerical value of the constant a in Eq. ( |5.3| ) comes from onr estimate 

aaf ( 247 / 9 V 3 - 23/2 -h 14/37r) 


of leading (~ /i^) perturbative mass correction. The approximation used in this estimate is similar 
to that proposed in . We will present this calculation elsewhere. 
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The mesons described above are excitations over the stable vacuum of the system 
(which is unique for h ^ 0). If m > 0 and h is sufficiently small, the system exhibits 
also an unstable “false vacuum”, a global resonance state whose (complex) energy is an 
intensive quantity, i.e. 

^meta = RFmeta{m,h) , (5.4) 


where R ^ oo is the spatial size of the system. The corresponding energy density Fmeta is 
a complex-valued quantity, and its imaginary part is interpreted as the decay probability 
density (the probability per unit volume and unit time) of the false vacuum. According 
to standard arguments (see [PU| , ^ , |4T[1 ) the resonance energy density h) coincides 

with the analytic continuation of the vacuum energy density F(m, h) in h, from positive 
to negative values of h, i.e. it can be written as 

.2 


m 


Fmetaim, h) = ^ log \h\ 15 ^metaiv) , U = m/\h\ 15 , 

Stt 


(5.5) 


where ^meta{y) is the function dehned in Sect. 3. According to (|3.18|) , at small h (i.e. at 
large y) the imaginary part 

r(m, h) = \h\^ Dmetaiy) (5.6) 

decays exponentially, see Eq. (|3.49|) . In fact, it was argued in that the asymp¬ 
totic (|3.18|) itself holds with exponential accuracy provided one replaces 2 ah AF{m, h) 
and m —>■ mq{m, h), where 


AF(m, h) = 3fJe Fmeta{rn, h) — F{m, h) ~ 2ah -f 0{h^) , (5.7) 


is the “string tension”, and mq{m,h) is the “quark mass”. Although it is not clear how 
to give a precise dehnition to the last quantity, it is natural to assume that at small h it 
admits an expansion (perhaps, an asymptotic one) in powers of h?, 

mq = m (d -f O j , (5.8) 

where q 2 is some constant, whose evaluation constitutes a separate and somewhat involved 
problem; at the moment we have only a preliminary numerical estimate, q 2 ~ 0.14(1) i. As 
the result, the expected asymptotic form of the imaginary part of the function 4>meta(?/) is 

Dmetaiy) ^V{y) (yQ + Viy~^ ^ -j y^Foo, (5.9) 

^ It results from estimating the three-quark contribution to the quark self-energy; we plan to 
address this problem separately. 
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where V{y) is defined in ( p.50|) , and 


— —2 

Vo = A = 0.216106..., Vi =-^^ + —I-^-0.14. (5.10) 

27r 2 4^1 

In TFFSA we actnally stndy the system in finite-size geometry, with the spatial com- 
pactification length R. In this sitnation the above stable mesons Mi correspond to a series 
of finite-size energy levels Ei{R) which asymptotically {as R ^ oo) behave as 

E,{R)=FR + M, + o{e-'^^^^y (5.11) 

where E = E{m,h) is the vacnnm energy density. Of conrse, there are no resonance 
states at finite R. If h is snfficiently small, the “false vacnnm” resonance of the R = oo 
system is very narrow, and at finite R it manifest itself as a pecnliar pattern of “near¬ 
intersections” of the levels Ei{R), see Fig. 6a in Sect. 6. The near-intersection of the levels 
Ei{R) and EiJ^i{R) occurs approximately at Ri^i+i = Mi/AF, where AF{m,h) is the 
same as in (13), and between the two near-intersections the level Ei{R) approximately 
follows the “resonance” linear law 


Ei{R) K, g?e Fmeta{m, h) R for Ri-i,i < R< R 


*,*- 1-1 


(5.12) 


This pattern (which was previously observed in [^, in the the TCSA study of IFT (|1.1|) ), 
is clearly visible in Fig. 6a in Sect. 6. Calculation presented in Appendix B shows that the 
separation between Ei{R) and EiJ^i{R) close to their mutual near-intersection follows the 
square-root law 


Eij^i{R) - Ei{R) K.mt^ \ {R - Ri,i+i)‘^ -F 


4TR 


mt‘^ a'- 


/ ’ 


(5.13) 


where F = F(m, h) is the resonance “specific width”, t = (2s^)^/^, and the constants a'i are 
the derivatives of the function cr(e), Eq. ( p.l6|) , taken at its zeroes e^. According to ( p.l7| ) 
and ([B.19|) , these constants have the small-^ expansions 


Ai'(-z,) / 2 19 ,4 , ^ 

BiFyr I “ 1400 * ) 


(5.14) 


Eq. ( |5.14|) will be used in Sect. 7.3 for the numerical determination of F(m, h) at 
small h. 
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6. Numerical Analysis 


In the numerical analysis of the Hamiltonian (|2.3D we used the technicaly simplest 
(albeit may be not optimal) truncation scheme based on the notion of level. We dehne 
the level of the state (|2.5a|) (or (|2.55|) ) as a half of the sum of absolute values of the 
half-integers ki (or the integers ni). At m = 0 this dehnition coincides with the standard 
notion of the level in TCSA. The truncation level L is the maximal level of states admitted 
into the truncated space. For the purposes of this paper we restricted our attention to 
the zero-momentum sector, i.e. the states with 'Yhi 'Thi equal to zero. For a given 
truncation level, the Hamiltonian (|2.3|) was numericaly diagonalized in this sector. We 
used different truncation levels ranging from 10 to 12; the dimensionality of the truncated 
spaces then ranges from 487 for L = 10 to 1186 for L = 12. This way some number of 
lowest energy levels Ei{R) was obtained for various real values of h, corresponding to rj in 
the interval [—5 : 5], as well as for pure imaginary h, corresponding to y in (|3.37|) ranging 
in the same interval. 


6.1. Finite-size levels. Qualitative picture 

The qualitative pattern of the resulting hnite-size levels Ei{R) is not very sensitive to 
the truncation level (except for some close proximity of the Yang-Lee point, see below), as 
long as R is not too large. Of course, higher truncation levels allow for better accuracy in 
quantitative estimates of its characteristics (energy density and masses). Also, as is usual 
with truncated Hamiltonians, the truncation effects become more prominent for larger 
values of R. In quantitative analysis we mostly used the data with R\h\'k <6. 


6.1.1. Real h 


The patterns of the lowest hnite-size energy levels Ei{R) for some real rj are shown 
in Figs. 6a-6c. In the range of R shown in these Figures the difference between the plots 
for L = 10 and L = 12 would not be visible. Qualitatively, the patterns in Figs. 6a-6c are 
the same as were observed in in the TCSA study of (|1 . 1|) , where their interpretations 
in terms of the particle spectra are discussed. In this case, we have just a few words to 
add. 


In full accordance with the McCoy and Wu scenario , for sufficiently large positive ry 
(i.e. for T < Tf, and small h) the held theory (|1 . Ij ) contains a large number of stable 
particles (“mesons”). At hnite R, these particles correspond to the energy levels Ei{R) 
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T1 =2.2 



ri =0.5 



R 

Fig. 66 


ri =-3.0 



Fig. 6: Plots of several lowest energy levels 
Ei{R), i = 0,1, 2,..., at different rj. 

a) rj = 2.2. Magnification in the corner shows 
the “near-intersection” of the levels Ei (R) and 
E2{R). 


b) 7] = 0.5. Three lowest excited levels Ei{R), 
E 2 {R), Es{R) are seen to approach straight 
lines ( ^.11 ) exponentially; these levels cor¬ 
respond to three stable mesons Mi, M 2 , Ms- 
The rest of the levels correspond to scattering 
states. 


c) rj = —3.0. Only one stable particle remains. 


which approach, at sufficiently large R, the straight lines F{m, h) R + Mi, where Mi are 
the masses of the mesons. This pattern of levels is clearly visible in the plot of Fig. 6a 
representing the case of ry = 2.2. For sufficiently large 77 the meson masses Mi are in 
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good agreement with (|5.2|) (see Fig. 5). The plot in Fig. 6a also shows a characteristic 
pattern of “near intersections” of the meson levels with the straight line 3?e Fmeta R, which 
is the hnite-size manifestation of the “false vacunm” already mentioned in Sect. 5. As is 
discnssed there and in Appendix B, the openings at these near-intersections are related 
to the imaginary part (^.6|) of the “false vacnnm” specihc energy, see Eq. ( ^.13|) . We will 
nse this relation in Sect. 8 as a way to estimate the resonance width at snfhciently large 
positive Tj. 

As rj decreases, the mesons with higher masses gradnally disappear into the 
“continuum” above the the stability threshold 2 Mi, where Mi is the mass of the lightest 
particle (see Fig. 66, where only three stable particle levels are apparent). At the same 
time the resonance “level” becomes fussy, in agreement with the expectation that its width 
parameter F grows. The process of depletion of the “meson” spectrum continues as rj de¬ 
creases to negative values. The third lightest particle disapears at 77 ~ —0.14, and the 
second one leaves the spectrum of stable particles at 77 = 772 , 

772 = -2.09(4). (6.1) 

The way these particles disapear is an interesting question which we hope to address 
elsewhere. For 77 below 772 only one particle with mass Mi is left, as is manifest in Fig. 6 c. 
The 77 -dependence of the hrst three masses was shown in Fig. 5. We postpone a detailed 
quantitative discussion of the particle spectrum to future publication. 

6.1.2. Pure imaginary h 

Much different patterns of the hnite-size energy levels appear if one takes h pure 
imaginary, which corresponds to 77 = y with real y. Although in this case the Hamil¬ 
tonian (|2.3|) is not an Hermitian operator, it has an explicit symmetry 

SHS = H^, ( 6 . 2 ) 

where the operator S just hips the signs of all states in the R sector in ( |2.5fc|) . As the 
result, the eigenvalues of this Hamiltonian are either real or come in complex conjugated 
pairs. 

For sufficiently large negative y the ground state energy Eq{R) of the truncated 
Hamiltonian remains real in a wide range of R, including the region where the linear 
behaviour ( |2.4|) is already clearly visible. This is also true for the hrst excited level Ei{R), 
which quickly (i.e. exponentialy) approaches the straight line E R + Mi, i.e. it behaves as 
a one-particle state with a real mass Mi. The situation is exemplihed in Fig. 7a, where 
the real parts of the hrst few energy levels are plotted. 
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Re E(R) / I h I , Re E(R) ^ ' h I 


y = -3.8 



y = -2.4 



y = -2.5 



R 

Fig. 7b 


Fig. 7 : Plots of several lowest energy levels 
Ei{R), i = 0,1,2 ,at pure imaginary h, i.e. 
at r] = y exp(47ri/15) with real y. When a level 
becomes complex, only its real part is shown. 

a) y = —3.8. The three lowest levels Eq(R), 
Ei{R), E 2 {R) are real for all R\h\'i^ G [0;6]. 
The level Ei{R) approaches the straight line 
E R + Mi exponentially; it corresponds to one- 
particle state. 

b) y = —2.5 (close to the Yang-Lee point 
Yq ~ —2.43). The three lowest levels remain 
real for all i?|h|T5 G [0:6], but the gap be¬ 
tween Eq{R) and Ei{R) narrows. 

c) y = —2.4 (on the other side of the YL 

point). The “nose” appears: Eq{R) and Ei{R) 
collide at R = R^ose ^ and be¬ 

come a pair of complex-conjugated levels at 


As y increases, the gap Mi becomes smaller (see Fig. 7b), and for y above —2.43 
a collision of the hrst two levels at some hnite R = Rnose{y) (the “nose”) is clearly 
visible; at i? > Rnoseiy) the levels Eq{R) and Ei{R) form a complex-conjngated pair, 
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as is seen in Fig. 7c. When y fnrher increases, the “nose” qnickly gets shorter, and when 
it is short enongh (as in Figs. 8a-86) it becomes evident that for large R 3> Rnose the 
complex conjngate levels Eq{R) and Ei{R) develop into the “straight lines” with com¬ 
plex slopes, Eq ^ E R, El (R) E* R. We interpret this as the vacnnm state acqniring 
a complex energy density E (obvionsly, in this case there are two “complex conjngate” 
vacna). Moreover, the next lowest levels E 2 {R) and Es{R) eventnaly collide and they too 
become a complex conjngate pair, both then approaching similar straight lines, shifted by 
some complex “mass” Mi, i.e. E 2 {R) E R + Mi , Es{R) E* R + Mj^; these levels 
are natnraly interpreted in terms of the complex-mass particle-like excitations over the 
corresponding vacna. 


y = 2.0 y = 2.0 




Fig. 8: The lowest energy levels Ei{R), i = 0,1,2,..., at y = 2.0 X a) Real 

parts of Ei{R); b) Imaginary parts of Ei{R). Two lower levels, Eq, Ei, collide at 
R = Rnose ~ 0.389 becoming at larger R a complex-conjugated pair. They quickly approach 
straight lines, whose slopes are interpreted as F and F*. Similar collision of the next two 
levels, E 2 , Es, occurs at i? ~ 2.28; they also develop into two complex-conjugated straight 
lines, with the same slopes as F and F*. The gap E 2 {R) — Eq{R) approaches a complex 
constant Mi, interpreted as a complex mass (correspondingly, the gap E^{R) — Eq{R) ap¬ 
proaches Ml). 
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As y becomes larger the number of such complex-mass particles increases. Higher energy 
levels form a rather complicated pattern. Their interpretation in terms of particle scatter¬ 
ing states (and even legitimacy of such interpretation) in general remains an open problem 
which we hope to address in the future. 

Appearence of complex F for y above the Yang-Lee point —Yq is well expected in view 
of the discussion in Sect. 3; F and F* are just the values of F(m, h) on two edges of the 
branch cuts in Fig. 1. We can say that the patterns in Fig. 7b and in Fig. 7c represent 
the situations at the opposite sides of the Yang-Lee edge singularity, which is located 
some place in between. Unfortunately (but not surprisingly), in the close proximity of this 
critical point the results from the truncated Hamiltonians appear to be more sensitive to 
the truncation level then elsewhere, especialy for larger values of R. In Fig. 9 we show the 
plots of the function Co{R) = ^ {Eq{R) — FqR), with 

Fo ^ 0.54743 (6.3) 

computed for different truncation levels, at y = —2.4295 (which we believe to be very close 
to the Yang-Lee point — Yq). As the truncation level increases it seems to develop the 
expected large-i? behaviour (|3.9|) . 


-0.39 


^ -0.4 

0? -0.41 
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^ -0.42 

Di 
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-0.43 


0.5 


y = -2.4295 


1 1.5 


.. , , , 

L= 10 
L= 11 
L= 12 
from Eq. (3.9) 


2 2.5 3 

R I h I 


3.5 


4.5 


Fig. 9: The function Co{R) = 6R[Eo{R) — 0.54743|/i|t5 A] / tt computed through 
the TFFSA using different truncation levels L = 10,11,12, at y — —2.4295xe^’^®A5, 
The dotted line is the expected asymptotic value —Ceff = —0.4. The dashed line 


represents the two leading corrections written in Eq. (^) with ag = —4.2. 
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6.2. Quantitative analysis 

The free energy F{m, h) was determined by measnring the slope of the R dependence 
of the gronnd sate level Eo{R) at snfhciently large R. There are two obstacles (well known 
in TCSA) to a precise determination of this slope. Theoretically, the linear behaviour (|2.4| ) 
becomes exact for large R 3> while at hnite values of R hnite-size corrections may 

be important. Although the leading hnite-size corrections are related to the masses Mi of 
the stable particles by well-known formula (see e.g. 0 ) 

Eo{R) = RF{m,h) -J" — Ki{MiR) + O (e-^ , (6.4) 

i 

where Ki{x) is modihed Bessel function, the subleading corrections depend on the (ge¬ 
nerally unknown) scattering amplitudes, and in most cases it is not clear how to take them 
into account. This makes it favorable to use larger values of R, where these subleading 
hnite-size corrections are less signihcant. However, in practical calculations the ground- 
state energy deviates from the linear behaviour (|2.4|) at large R as the result of the state 
space truncation (the “truncation ehects”). Although this deviations are too small to be 
visible in Figs. 6 - 8 , they substantially reduce the precision of quantitative analysis. The 
truncation ehects bring an upper limit to the values of R at which precise quantitative 
analysis can be made. 

To reduce the truncation ehects, the ground-state energy functions E^\r) computed 
with the truncation levels L = 10,11,12 was extrapolated to L — oo. The extrapolation 
was made by htting the formula 

Ei^\R) = so{R) + Si (A) (6.5) 

to the computed L dependence (at hxed i?), via so,si,S 2 , and accepting so{R) as the 
extrapolated Eq{R). Although this formula does not have theoretical justihcation, the ex¬ 
trapolation procedure appears to work well for all real as well as imaginary h, corresponding 
to 1 77 1 < 5, substancially improving the linear shape of Eo{R) (e.g. for R 20M^^ the 
deviation from the linear behaviour is reduced by at least one order of magnitude). 

The extrapolated function Eo{R) was then subjected to a hnite-size analysis. The 
masses Mi of the stable particles were estimated from the higher levels Ei{R), and the 
expected leading hnite-size corrections (the sum in (| 6 ^)) were subtracted from Eo{R). 

For negative ry in Table 4 we subtracted all one-particle corrections and then used 
the function foR + fi exp(— 2 / 2 A) to ht (around R ~ 4M^^) the remaining part via 
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fo: fi: f 2 , with the best-ht value of /o accepted as F{m, h). For positive r] we subtracted 
instead the corrections from the two lightest particles and then used qqR + gi Ki{g 2 R) 
to ht the resulting curve (around R ~ 7 M^^). The overall quality of this analysis was 
controlled by observing the proximity of the best-ht value of /2 and (72 to the expected 
values 2 Ml and M 3 , respectively (the deviation varied with 77 , but never exceeded 20%). 

The above procedure of the hnite-size analysis is rather labor consuming. In order 
to obtain estimates of the coefficients we used similar procedure in the small inter¬ 
val [—0.64 : 0.64] around zero but we only subtracted the lighest particle correction before 


-1 


htting to go R + gi Ki (772 R) around i? ~ 8 . 

We believe that the data for F(m, h) obtained through this procedure are exact to six 
signihcant digits. 

Finally, massive data used in Fourier analysis in Sect. 8.2.1 was obtained using some¬ 
what simplihed procedure, in which only the hrst term in the sum ( |6.4|) was taken into 
account, i.e. no subtraction was performed. Specihcally, the function goR + gi Ki{g 2 R) 
was htted to the extrapolated data Eo{R). This procedure produces the estimates which 
agree with the above more elaborate analysis to hve signihcant digits. 

Similar simplihed analysis was used in the case of pure imaginary h. In this case, 
for y > —To, complex values of the htting parameters go, gi, 772 were admitted. Exception¬ 
ally, for the values on the side y < —Yo presented in Table 5, we subtracted the leading 
hnite size correction and then htted to goR + gi exp(— 772 /?). Special was the vicinity 
of the critical point —Fq. With y approaching this point the values of 772 decreased and 
overall quality of the ht deteriorated. Of course this is not surprising since near the critical 
point the correlation length M^”^ is diverging, pushing the domain of validity of the expan¬ 
sion ( |6.4| ) to i? —> cxD, where our analysis would be severely spoiled by truncation ehects. 
Nonetheless, it was possible to obtain reasonably accurate data (at least four digits) for 
^imh{y) with y < —2.43. On the other side of the critical point, for y > —To; the situation 
is complicated by the rapid growth of the “nose”. In this region we could obtain similar 
accuracy only for —2.1 < 77 , and our data for —2.3 < y < —2.1 may be less accurate. In 


close vicinity of the critical point a hnite-size analysis based on the ehective action (|3.7|) 
would be more appropriate, and we are planning to apply it in the future. 
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7. Numerical results 


7.1. Scaling function. Real h 


Our results for $(? 7 ) at real values of ry are shown in Fig. 10 (see also Table 4 for some 
numerical values). As expected (see Sect. 3), the function appears to be perfectly smooth 
everywhere exept for the point 77 = 0 , where the singularity of the form — 77 ^ logTy^/STr can 
be observed. The hrst few coefficients of the 77 -expansion ( |3.29[ ) can be readily extracted 
from this data. After subtracting the above singular term, we have htted polynomials 
of various degree to the data in the interval [—0.64 : 0.64]. The ht appears stable for 
the hrst eight coefficients, $0 to $7. The resulting estimates of these coefficients are 
displayed in Table 3. The estimates for the coefficients $0 and $1 are very close to the 
predictions (|3.31| ) and ( p.33|) ; this result for $1 can be viewed as numerical verihcation of 
the exact expectation value (e)|T-=o obtained in [J3|. It turns out that the hrst few terms 
of the series (|3.29|) provide rather good approximation for $( 77 ) with I 77 I < 2, as shown 
in Fig. 10. Comparison with the known terms of the expansions (|3.35 ) and ( 3.36 ) is also 
presented in that Figure. 
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Fig. 10: The scaling function ^{rj) for real 77 (solid line). Dashed lines show how 
this function is approximated by few leading terms of the expansions (3.29), (3.35) 
and ( 3.36 ). 
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Fig. 11: The TFFSA data for the function ^imh{y)i Eq. ( p.37|) (empty and full 
bullets). The imaginary part vanishes for y < —Yq. Various dashed 

lines show how this function is approximated by its expansions ( 3. 291) , (|3.35D and 
( 3.36 ). 16 and 17 are the first seven (including $6 2/^) and eight (including ^jy'^) 
terms of (H). Corresponding expansions of the real part are the lines R6 and 
R7. One, three, and four terms of the asymptotic expansion (|3.36 ) are shown as 
LTl, LT3, and LT4, respectively (obviously, only even (odd) terms contribute to 
the real (imaginary) part). Finally, HT12 is the plot of six leading terms (including 
Gi 2 term) of the low-T expansion ( ^.35| ). In all cases the coefficients are taken (or 
computed) from Tables 1, 2 and 3. 


1.2. Scaling function. Imaginary h 

Our results for ^imhiy) are shown in Fig. 11, where a singularity aX y — —Fq: 


Vo ~ 2.4295, 


(7.1) 


is clearly visible. The function ^imhiy) is real for y < Yq, and it develops an imaginary 
part for y > Yq. A comparison with the hrst few terms of the expansions (|3.35|) and (|3.36| ) 
is also made in the Fig. 11. We observe that the hrst six terms of (|3.35|) approximate 
^imh{y) very well for y < —2.7. This is not surprising since, according to our discussion in 
Sect. 3, the expansion (3.35|) is expected to converge for all y <Yq. In this Figure we also 
compare the data with the hrst eight terms of the expansion (p.43|) (and similar expansion 
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for the real part of ^imhiv))- One can notice that these expansions provide a remarkably 
good approximation of ^imhiv) for \y\ < 2. We take this as an indication that it is the 
Yang-Lee singnlarity at ?7 = —Yoe^^ which determines the domains of convergence of 
the small ry expansion (|3.29|); in other words, it snpports our expectation that the function 


$(? 7 ) is analytic in the whole disk I 77 I < Yq- 

Our data agree with the expected form (|3.40| ) of the singularity ai y = —Yq. The 
imaginary part of Eq. ( |3.39|) , receives no contribution from the regular part A[y) 

in (|3.40| ), and therefore it should exhibit the singular behaviour (|3.4CI|) most prominently. 
The log-log plot of this function in Fig. 12a clearly confir m s the value 5/6 of the exponent 
in ( b.40|) predicted in The ratio Am^imh{y)/{y + Yq)^ is shown in Fig. 126. 
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Fig. 12: TFFSA data for Amh(y), Eq. (13.391) . a) Log-log plot of Amh(y) (bullets). 
The dashed line has slope 5/6; 6) Ratio Di^h{y)/{y -\- Yq)^^^ (bullets), and the graph of 

(dashed cnrve). 


Amh(y)/(y + Yo)®/®, with Amh(y) from Eq, 


This plot reveals the “hne structure” of this function near the Yang-Lee branching 
point, which should be attributed to the subleading singular terms in (|3.4C1|) . Our data is 
not very precise when taken close to the critical point — Yq (see Sect. 6 ), and allow for only 
rather rough estimate of the leading coefficients in the expansion (|3.40|) , (3.41), 


5o =-1.37(2), 
Ao = 0.5474(3), 


Bi = -0.75(5), 
Ai = 1.06(4), 
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Co = 0.38(7), 
A 2 = 0.3(1). 


(7.2a) 

(7.26) 



























These estimates are obtained by fitting a few leading terms of (3.41) to the data in varions 
intervals of y between —2.2 and —1.5; the error reflects the dependence on the htting 
interval as well as on the nnmber of the terms nsed. From (|7.2a|) we can obtain an estimate 
of the most important parameters in ( p.8a|) , ( p.8fi|) , 


Ai = 3.10(6), CKO = -4.9(9). 


(7.3) 


Independent (and somewhat better) estimate of the parameter ckq can be obtained di¬ 
rectly from the hnite-size ground-state energy Eo{R) at the Yang-Lee point y — — Yq- 
By htting the expected form (iD to the data in Fig. 9 in various regions between 
R = 1.7 |h|“^ and R = 2.8 |h|“A we found 


CKO = -4.4(4), (7.4) 

(where again the error estimate is from the dependence on the htting region). The quality 
of the ht is also shown in Fig. 9. 


7.3. False vacuum resonance width 


As was discussed in Sect. 5, the hnite-size behaviour of the meson levels close to their 
“near-intersection” points is related through (|5.13|) to the width F of the false vacuum 
resonance, and can be used for numerical determination of Dmeta{y)- Presently, there 
are two limitations to this approach. First, at small h the “near- intersections” occur at 
large R, where truncation ehects are more signihcant. In our calculations, a satisfactory 
stabilization of the diherence E 2 {R) — Ei{R) at L = 12 was achieved only for y < 3.5. For 
larger rj the situation can be somewhat improved by adding a number of two-quark states 
of higher levels (see below). On the other side of large h, a limitation comes from our 
unsufhcient knowledge about the coefficients a' in (|5.13|) ; the best we have at the moment 
is the hrst three terms (|5.14D of their small-h expansions. Nevertheless, one can note that 
these three terms were obtained in Appendix B within the same approximation as the 
three correction terms in the large-expansion of the meson masses, Eq. ( b.2|) . Good 
agreement of the last expansion with the data at ry > 1.7 makes us believe that reasonably 
accurate estimate of the function D.^netaiv) can be obtained within this expansion in the 
same region. 

We have performed this analysis in the range 1 < ry < 4.5 using the separation 
E 2 {R) — Ei{R). The eigenvalues Ei{R), E 2 {R) were computed in “two-quark extended” 
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truncated space (2.5) with L = 12. Namely, besides all states in (2.5) with L < 12, the 
two-quark states \ 0)ns and \ 0)r with the levels 12 < L < 100 were also 

admitted. While such extention does not bring any appreciable difference for rj < 3.5, 
it does snbstantially improve the resnlts for 3.5 < rj < 4.5. For all valnes of rj in the 
range 1 < ry < 4.5 the separation E 2 {R) — Ei{R) was in excellent agreement with the 
expected square-root behaviour ( b.l3|) , and it allowed ns to estimate the values of the 
width parameter F. The results for the ratio Dmetaiv)/V(y) (where Dmetaiv) is related to 
F through (|5.6|) , and V{y) is the fnnction ( |3.50|) ) are shown in Fig. 13. The ratio shows a 
clear tendency to approach at large y the expected valne Vo = 0.2161.... Moreover, when 
all three terms in (|5.14|) are taken into acconnt, for y > 2.3 it agrees well with the first 
correction term in ( b.9|) . By direct fitting the expansion ( |5.9|) to this data we have obtained 
the following estimates of the coefficients there, 

W = 0.2161(4), Vi = -0.136(4). (7.5) 



Fig. 13: Data for the ratio D^eta{y)/V{y) obtained through the false vacuum 
resonance width analysis (Sect. 7.3), using Eq. ( b.l3| ) with three leading terms in 
Eq. (5.14). The data obtained using “two-quark extended” truncation space with 
the truncation level L = 12 is shown as empty bullets (0). The data obtained 
using two-quark states only is shown as full bullets (•). The solid line is the plot 
of Eq. (|8.11D and the dotted line is the value Vq = 0.2161. 
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In [|2] an attempt was made to determine the h ^ 0 behavionr of F using two-quark 
approximation for the resonance wave function. While correctly reproducing the exponent 
in (|3.18|) , this calculation yields somewhat greater pre-exponential factor rrah/18 (which 
would lead to Vq = tts/IS ~ 0.237 in our notations). As we explain in Appendix B, 
the two-quark approximation breaks down when the quarks are deep inside the classically 
unaccessible region, and therefore it is not justihed in the resonance wave function calcu¬ 
lation even at h ^ 0, even though in this limit the resonance decay goes predominantly 
through the two-quark channelS. To verify this being the reason for the above discrepancy, 
we have repeated the calculation described in this Section, this time with the eigenvalues 
Ei{R), E 2 {R) computed purely within the two-quark approximation. Specihcally, in com¬ 
puting these eigenvalues, only two-quark states with the levels < 100 were admitted. As 
is seen in Fig. 13, the ratio Dmetaiv)/Viy) obtained this way is hardly consistent with the 
correct y ^ oo asymptotic value 0.2161.... 


8. Analyticity 


8.1. High-T domain (HTW in Fig. 4) 

In principle, our approach allows for direct computation of the function $(?7) at com¬ 
plex values of r] in the wedge HTW in Fig. 4. However, this data would be redundant 
because in view of the analytic properties of the scaling function discussed in Sect. 3, the 
values of ^{rj) can be computed through the dispersion relation ( p.3|) once the function 
Dimh{y) = is known. This dispersion relation applies to the function ^{rj) 

with ry in the wedge HTW in Fig. 4., where it can be rewritten as 


$(? 7 ) 


15 Dimh{-y)y* dy 

4 7r io y^ + 


( 8 . 1 ) 


Since accurate direct data on Dimh{y) in the vicinity of the Yang-Lee point (—Tq < y < 2.1) 
is not available at this moment, some extrapolation is needed in order to perform the 


(numerical) integration in ( p.l|) . We have tried various extrapolations, and the results 


^ To avoid misunderstanding let us stress here that in Appendix B we apply the two-quark 
approximation to the wave function in the “meson domain”, not deep under the barrier. Instead 
of solving the problem in the last region, we take a phenomenological approach based on the 


Breit-Wigner formula (B.23). 
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of integration in (8J.) usnaly agree with the our direct data on $(? 7 ) to at least three 
signihcant digits. The best results were obtained with the following approximation 


Dimh{y') 


-Y„<y< -1.57; 
DLyiv) for -1.57<i/<0. 


(8.2) 


where 


with 


OM(y) = -lBo(v + Yo)i + ^Bi{v + Yo)* + ^C„{y + Y„)i , 


Bo = -1.3693 , Bi = -0.74378 , Cq = 0.42446 , 


(8.3) 


(8.4) 


and D^^{y) is just the hrst eight terms of the y-expansion, 


D^mhiy) = $0 + y + ■ ■ ■ + $7 y ■ 


(8.5) 


with related to 4*^ from the Table 3 as in (|3.44|) . Here we accept ( [7.1|) as the value of Yq, 
while the coefficients Bo,Bi,Co in ( ^I^) are obtained by exact matching (|8.3|) with (|8T| ) 
at y — —1.57 (up to the hrst derivative), and hne-tuning the remaining one parameter 
to achieve better agreement with the data. Note that these coefficients agree with the 

should not be considered as 


estimate ( [7.2a|) within the stated accuracy. However, 
better estimates of the coefficients of the expansions (3.416)-(3.41c), since (|8.3| ) ignores 
higher subleading terms in (|3.40|) . The quality of the approximation (|8.2|) is demonstrated 
in the Tables 4 and 5, where the results from (|8.1|) are compared to our data, both for 
real 77 < 0 and y = ye^ ,y < —Yq. We also presented in Table 1 the results for the 
hrst few coefficients G 2 n in ( P.l|) obtained through the dispersion relation (^) with the 
use of o. The agreement with the previously known estimates is reasonable, baring in 
mind relatively low quality of the approximation (|8.3|) ; we believe it will improve once the 
structure of singular expansion (|3.40|) is understood in greater quantitative details. 
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$( 77 ) 

From 

TFFSA data 

From high-T 
disp. relation 

From low-T 

disp. relation 

From extended 
disp. relation 

*(-5) 

-0.1092101 

-0.1092092 

— 

-0.1088626 

<l>(-4) 

-0.1592682 

-0.1592643 

— 

-0.1589421 

<I>(-3) 

-0.2529928 

-0.2529887 

— 

-0.2527417 

*(- 2 ) 

-0.4413450 

-0.4413249 

— 

-0.4412136 

<I>(-1) 

-0.7839665 

-0.7839668 

— 

-0.7839576 

<l>( 0 ) 

-1.1977330 

— 

— 

-1.1977320 

<I>( 1 ) 

-1.3898410 

— 

-1.3898417 

-1.3898063 

4(2) 

-1.4930558 

— 

-1.4930566 

-1.4929849 

<l>(3) 

-1.5642732 

— 

-1.5642736 

-1.5641727 

4.(4) 

-1.6188506 

— 

-1.6188510 

-1.6187275 

4(5) 

-1.6632483 

— 

-1.6632485 

-1.6631076 


Table 4: Numerical values of 4*(r?) at selected real points. The first column shows 
direct data obtained through TFFSA, as explained in Sect. 6.2; we believe these 
nnmbers are exact to six significant digits or better. The resnlts from high-T and 

with the use of approximations 


low-T dispersion relations, Eqs. (| 8 .lD and 
(8.2) and ( ^.14 ), are presented in the second and the third columns, respectively. 
The last colnmn contains results from extended dispersion relation, Eq. (4.8), with 
approximation (8 .21 ) used. 


^imh{ 4.5) 

4’im?i(~3.5) 

$zm/i(-3.0) 

^imh{ 2.5) 
$^^,,(- 2 . 4295 ) 


From 

TFFSA data 

0.1116003 

0.1349420 

0.1674520 

0.2155503 

0.2946233 

0.4729294 


From high-T 
disp. relation 

0.1115969 

0.1349380 

0.1674469 

0.2155426 

0.2946092 

0.4728305 

0.5475373 


Table 5: Numerical values of ^imhiv), Eq. ( 3.37 ), at some real y < —Yq. Direct 
data from TEESA (good to four digits or better, see Sect. 6.2) are in first column, 
and results from high-T dispersion relation with (8.2) in second column. 
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8.2. Low-T domain (LTW in Fig. 4) 


8.2.1 Fourier analysis 

Numerical analytic continuation into this domain was performed using Fourier 
analysis in the logarithmic variable log 77 . According to (|3.29| ) and ( |3.36|) the function 
ryif ^$( 77 ) — Gi 77 ^^ decays sufficiently fast (exponentially in the above logarithmic vari¬ 
able) both at ?7 —> 0 and at 77 —> cx). Its logarithmic Fourier transform 

^(cj) = r 77 ^-^"^ ($( 77 ) - Gi 775 ) ^ (8.6) 

Jo V 

was evaluated numerically from our data for $( 77 ). In principle, using this f un ction the 
values of ^lowiv), Eq. (B, along the rays 77 = with real y G ( 0 ,-t-cx)), can be 

computed as the inverse Fourier transforms 

/ OO 7 

(8.7) 

-00 


The function ^lowiv) is analytic along the ray provided the integrand in (|8.7|) decays 
faster then any power of w. The standard analyticity assumption for the wedge LTW 
in Fig. 4 is equivalent to the statement that |T(c<;)| < exp (—^ 7 r|c<;|) for |ct;| —>• cx. Of 
course it is not possible to achieve a rigorous conclusion about true large to asymptotic of 
this function using our limited-precision data. For uj > Umax ~ 5.5 the numerical Fourier 
transform T(ci;) is dominated by a small (~ 10 “®) random-looking function which should 
be interpreted as the Fourier transform of the errors associated with our numerical data 
for $. For this reason the numerical Fourier transform \l/(c<;) hardly carries any usefull 
information beyond the window |c<;| < Wmax- At the same time the function \l/(c<;) inside 
this window appears to be decaying signihcantly faster then the above exponential, which 
is consistent with the standard analyticity assumption. One can actually use ( |8.7|) to 
compute this function in LTW with reasonable accuracy. Obviously, in order to do that 
one has to introduce some cut-off to eliminate the high frequency part |ci;| > Wmax of the 
integral (| 8 . 8 |) dominated by the noise; it is expected that only hne details of the resulting 
function would depend on the specihcs of the cutoff procedure. As an example (which will 
turn out to be usefull in our analysis in Sect. 8.3) we show in Fig. 14 the shape of the 
function 4>zoi(;(h) along the ray 77 = y > 0, computed this way. Any changes induced 

by changing the cutoff procedure are by far too small (< 10 “^) to be visible in this plot. 
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y 

Fig. 14: The scaling function ^lowiv) along the axis rj = y e~ obtained 

through numerical analytic continuation (Sect. 8.2.1). Real and imaginary parts 

_ 2712 

of e 3 ^iow{yG 3 ) are shown. The y < 0 part of this axis lays inside the shadow 
domain SHD in Fig.4, close to the axis OA. The data here is not very accurate. Nev¬ 
ertheless, an anomaly around y — —2.5 (presumably associated with the proximity 
of the YL singular point) seems to be developing. 


More interestingly, one can apply the relation 


D 


meta 


(y) — — y 1® cosh(87rci;/15) \l/(ci;) 


du 


( 8 . 8 ) 


(which is equivalent to the dispersion relation (|3.15|) ) to compute the function (|3.46|) . As 
expected, the resulting inverse Fourier transform (|8.8[ ) shows rapid decay at large y (see 
Fig. 15a). However, instead of monotonous decay ( p.49|) , the cutoff integral (|8.7|) exhibits 
small oscillations whose shape is very sensitive to the details of the cutoff procedure. 
Therefore, we consider the results for the function D^etaiy) obtained this way only reliable 


for y < 2.0. To see how these results comply with the expected asymptotic form (^) we 


have plotted in Fig. 156 the ratio Dmetaiy)/^(y), where V{y) is the exponent (|3.50|) . 
The error bars in this plot represent the spread of the data obtained with different cutoff 
procedures. From (^T^), this ratio is expected to approach the constant 0.2161... as y —> cx). 
As is seen from the plot, the ratio reaches a minimum 0.153 at y ~ 1.15, and then shows 
a tendency to increase towards the expected value. 
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y 

Fig. 15a 



Fig. 15: a) The function Ilrneta(y)i Eq. (|3.46|) , obtained through numerical analytic con¬ 
tinuation in Sect. 8.2.1; b) The ratio Ilmeta(y)/E(y) (where V{y) is defined in Eq. ( 3.50D ). 
The data from numerical analytic continuation (Sect. 8.2.1) (shown as error bars) vs. the 
data from resonance width analysis (Sect. 7.3) (shown as empty bullets). The solid curve is 
the plot of Eq. ( p.l4| ). 


Unfortunately, as was explained above, for y > 1.9 the data obtained by this method 
becomes unreliable, and presently neither the true large y asymptotic nor the way it is 
approached can be established on their basis. But fortunately, the analysis of Sect. 7.3 
apparently provides a complimentary set of data. We replot in Fig. 156 the data from the 
resonance width analysis (identical to those in Fig. 13) to demonstrate its agreement with 
the data obtained by the Fourier analysis. 

In principle, one could try to use this numerical Fourier transform \l/(ci;) to continue 
further, into the shadow domain, by taking |q;| > Stt/IS in (|8.7|) . We found that while for 
|q:| < 37r/5 the intergal in (|8.7|) still receives dominating contribution from uj well within 
the window |c<;| < Umax (and hence is not too sensitive to the cutoff procedure), for larger 
a the domain around Umax becomes more signihcant. The result of such computation 
with a = —2Tij3 is shown in the left-hand side of Fig. 14; its accuracy is hardly better 
then 10“^, and may be worse. For |q;| > 27r/3 the results are too sensitive to the cutoff 
procedure to be taken seriously. In other words, at present the accuracy of our direct data 
about ^low iv) is not sufficient to support (or unvalidate) analyticity of this function in the 
whole of the shadow domain. In Sect. 8.3 below we use another approach, based on the 
extended dispersion relation, Eq. ([4.8|), to verify the extended analyticity. 
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8.2.2 Low-T dispersion relation 


According to (|3.15| ), in the low-T wedge —Stt/IS < arg(? 7 ) < Stt/IS (the domain LTW 
in Fig. 4) the scaling fu n ction $( 77 ) is represented in terms of the fnnction Dmetaiv) as 


$(77) = Gi 778 


1577 ' 


8 TT 


Dmeta{y') V ® dy 

_ 15 _ 15 

y 8 + ?7 8 


(8.9) 


It is interesting to check how the above information abont Dmeta{y) agrees with our direct 
data on this scaling function. Inspection of the plot in Fig. 156 shows that the function 
D-metaiy) is Well approximated by the hrst eight terms of its power series expansion (|3.47|) , 
i.e. by the polynomial 


ddmetaiy) — Dq + Diy + - \- , 


( 8 . 10 ) 


whith the coefficients Dn computed from $n’s in Table 3, all the way up to y 
plot also shows that the above large -77 approximation, 


£>me.„(!/)= (0-2161-0.136!;- 


V{v), 


1.1. This 


( 8 . 11 ) 


seems to provide rather accurate approximation of this function for y > 2.5. The two 
curves corresponding to (|8.10| ) and (|8.11| ) intersect at 77 ~ 1.59. We have checked that 
even if one just takes (|8.1CI|) and (|8.11|) as the approximations for Dmetaiy) at 77 < 1.59 
and 77 > 1.59 respectively, the dispersion integral (|8.9|) reproduces correctly at least four 
signihcant digits when compared with our direct data at real positive 77 . Of course such 
accuracy is related to the fact at 77 ~ 1.59 the function Dmetaiy) is already rather small. 
In fact, much better approximation can be achieved if one uses some interpolation between 
( ^. 10 |) and (|8.11|) in the intermediate region. We have obtained very good results with the 
hfth-order polynomial interpolation 


DtSiv) = 


Po + Pi {y - vi) ^— + Paiv- yiY 


V{v), 


( 8 . 12 ) 


with 771 = 1.1 and the coefficients Pq, ... hxed by the conditions that ( 8 . 12 ) matches 
smoothly (up to the second derivative) the approximations ( | 8 . 1 C 1 |) and (| 8 . 11 j ) at 77 = 1.1 
and at 77 = 2.72, respectively. 


Po = 0.1532863 , Pi = 0.0002611, P 2 = 0.0499557 , 
P 3 = -0.0209331, P 4 = -0.004408 , P 5 = 0.002672 . 


(8.13) 
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The dispersion integral (|8.9|) with 


1 


for 

0 < 7 / < 1.1 


Dmetaijj^ ~ > 

\ ^interp/ 5 

for 

1.1 < 77 < 2.72 

(8.14) 

1 

1 

for 

2.72 < y < 00 



reprodnces onr direct data very accnrately (with six signihcant digits, i.e. essentially 
within the estimated accnracy of the direct data), as is demonstrated in Table 4. In this 
computation we have used the exact value (|3.12| ) of the coefficient Gi in (|8.9|) , and again 


the value (|7.1|) of Yq was assumed. 

The approximation ( 8.14 ) can be also used to estimate the coefficients of the 
assymptotic expansion ( 3.11 ). It follows from (|3.13|) that 


15 


G^ = (-Y+^ — 

^ ’ 87r 


_ (15n-24) 

D{y)y 5 dy 


for n > 1 


(8.15) 


Moreover, the coefficient Gi itself must coinside with the integral 

15 


Gl — — / {Dmeta 

StT /n 


iy) - Do)y ‘dy, 


(8.16) 


lest analyticity of $(? 7 ) at 77 = 0 be violated. With our approximation ( ^.14|) the integral 
( p.l6|) reproduces the coefficient Gi almost exactly (if fact, we have used this as a condi¬ 
tion in choosing the interpolation interval in (|8.14|) ), and few higher coefficients obtained 
through (|8.15|) are compaired with previous estimates in Table 2. 


8.3. Extended Dispersion Relation 


Equipped with the results in Sects. 8.1 and 8.2, we made a preliminary analysis of 
analyticity in the shadow domain in order to validate the extended analyticity conjecture 
formulated in Sect. 4. If one assumes the extended analyticity, the estimates in Sect. 8.1 
and 8.2 make it possible to actually reconstruct, with reasonable precision, the behaviour 
of the function $( 77 ) along the shadow side of the Yang-Lee branch cut in Fig. 4, thus 
hnding an approximation for the discontinuity fu n ction 


A{y) = i^imh{-y) - ie *15 


(8.17) 


which enters the extended dispersion relation (^4.8|) ; the expression ( p.l7|) is a simple con¬ 
sequence of (|4.3|) and (b:.l|) . The extended analyticity is then verihed by checking this 
dispersion relation against direct data. 
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Let us note that, on one hand, the results of Sect. 8.2 allow one to hnd the large-y 
asymptotic form of the function ^iow{y ^)- Indeed, the second (analytic) form of the 


Eq. (|3.46|) which dehnes Dmetaiv) also dehnes its analytic continuation to complex values 
of its argument. The extended analyticity is equivalent to the statement that Dmeta{z) is 
analytic in the wedge —tt/S < arg( 2 ) < tt/S. Under this assumption, it follows from (|3.46| ) 
that 




low [ye ‘ 15 = e 


' 15 ' 


^iow{ye3 )+2ie " Dmeta [y e 


^7^ 

5 


(8.18) 


For real positive y, this equation expresses the desired values of ^lowijli) on the shadow 
side of the branch cut in terms of its values well inside the domain LTW in Fig. 4, and 
also through the analytic continuation of the function Dmetaiy)- The hrst term is already 
under good analytic control, see Sect. 8.2.1. In fact, for sufficiently large y the function 
^iow{y e^) is approximated very well by hrst few terms of its asymptotic expansion (|3.36|) , 


as illustrated in Fig. 14. Specihcally, for y > 3.0 the hrst twelve terms of this expansion 
yield an accuracy better than ~ 10“'^. On the other hand, the large-y behaviour of the 
second term in ( p.l8| ) follows from the asymptotic law (|5.9| ). With this, Eq. (|8.18|) yields 
approximation for the second term in Eq. (|8.17|) , valid at sufficiently large y. This readily 
translates into corresponding large-y approximation for the function A(y), since the large-y 
behaviour of the hrst term in (|8.171) follows from (|3.385D and (|3.1|) . Finaly, one can check 
that for y > 3 the accuracy of at least ~ 10“^ can be achieved by taking the hrst six terms 
(i.e. including G 12 term) of the expansion (|3.1|). Therefore the following formula 


A(2/) 


12 

A (y) = 

k=l 


Gu 


e 8 Gi 


16-15fc 

y ® 


-\-2 [Vq -\-Vi e s y 8 1^8 exp 


_ 37r2 15 37rz I 

-^e ® ® r 


(8.19) 


is expected to approximate the discontinuity f un ction well for sufficiently large y. In ( 8.It ) 
y, stands for 7r/2s, and it is understood that G' 2 n+i = 0. 

On the other hand, when y is close to Yq the relation ( |4.2| ) is usefull. The approxima¬ 
tion (|8.3|) for Dimhiy) then suggests the following approximation 

My) ~ A<(y) = y‘^/4 

-0{y -Yo)(^l Bo (y - To) ^ ^ M e-M ^ (^^ _ y^) I j , 

( 8 . 20 ) 

which is expected to be valid at j/ ~ To- In (|8.20|) 9{x) is the step function. 
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Fig. 16: Plots of real and imaginary parts of A^(y), Eq. ( |8.20D (shown as solid 
lines), and A^(y), Eq. (8.19|) (dashed lines). 


The plots of (|8.19|) and (|8.20|) are shown in Fig. 16. In this plot the nnmerical val- 
nes ( |5.4|) of the coefficients in ( ^.201 ) are taken. Also, the coefficients Vq and Vi in (|8.iy| ) 
are identical to those in (|8.11|) , and the coefficients Gk and Gk are taken from the Tables 
1 and 2. In fact, contributions of the terms with A: > 5 in the sum in (|8.19| ), as well as the 
differences in the choice between the columns in that Tables, are too small to be visible 
in this plots. One can see that the approximations (^.19|) and ( ^.20|) match each other 
very well in the intermediate region of y around 3. This match itself can be viewed as a 
support to the extended analyticity conjecture. Even more convincing result is obtained 
by evaluating the extended dispersion relation (|4.8|) with the use of these approximations. 
We have evaluated the integral (|4.8|) with 


A(?/) = 



for 0 < y < 3 
for 3 < y 


( 8 . 21 ) 


for some real values of rj. In this calculation specihc numerical values of the coefficients 
G 2 n from the second column in Table I and Gk from the second column in Table 2, as well 
as the values ( b.II| ) of Vo,Vi, were used. The results for the function ^{rj) are compared to 
our direct data in Table 4. Also, the hrst coefficients 4*^ obtained through (4.10), within 
the same approximation, are presented in Table 3. Although the agreement is not very 
accurate (the values of 4’(ry) agree with the direct data with only four digits, which is well 
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below the estimated accuracy of the latter), we take it as a strong support of extended 
analyticity. The relatively low accuracy should probably be attributed to low quality of the 
approximation ( ^.201) ; we believe it would improve when better estimates of the coefficients 
in (3.416)-(3.41c) are obtained, and higher subleading singular terms in (|3.40|) are included. 


9. Discussion 


In this work we have used numerical data obtained by TFFSA to study the scal¬ 
ing f un ction $(ry) associated with the Ising model free energy in its critical domain. We 
estimated some parameters of this function, including location of the Yang-Lee singular¬ 
ity (|7.1|) and leading amplitudes (|7.2a|) , (|7.25|) of associated singular expansion (|3.4C1|) . We 
also examined the analyticity of the scaling function. Our analysis conhrmed the standard 
analyticity assumption (see Sect. 3), and we have found strong evidence in support of the 
“extended analyticity” conjectured in Sect. 4. Here we would like to explain why we regard 
proving (or invalidating) this extended analyticity as an important issue. 

As is known, in the classical theory of phase transitions, which does not incorporate 
fluctuational effects, the hanger’s singularity of the low-T free at ih = 0 is not present. 
Instead, analytic continuation of the classical low-T free energy shows a branch-cut sin¬ 
gularity at some hnite negative value oi H, H = known as the “spinodal point”. 

This singularity appears because at sufficiently large H {\H\ > Hqp) the metastable state 
becomes unstable against classical decay. When fluctuations are taken into account the 
free energy becomes singular (albeit weakly singular) at ih = 0, since the metastable phase 
is always prone to decay through nucleation. But what happens to the spinodal singu¬ 
larity? There are convincing arguments (see [^) to that in the presence of fluctuations 
the (analytic continuation of the) free energy cannot have singularities at real negative H. 
To our opinion, the most plausible possibility is that in fluctuational theory the spinodal 
singularity does not just disappear but instead gets pushed under the hanger’s branch 
cut, as is shown schematicaly in Fig. 2. In other words, the nearest singularity under the 
hanger’s branch cut could be associated with the “ideal limit of metastability” where 
a transition from a decay through nucleation to faster “classical” mechanism of decay takes 
place. The transition can be more or less expressed depending on how far under the cut 
the singularity is located. The above speculation is made under assumption that the ana¬ 
lytic continuation of the free energy is a quantity at all relevant to termodynamics of the 
metastable state, even away from the domain of vanishingly small H. It also ignores all 
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possible kinetic aspects of the problem, certainly important in realistic situations. Nonthe- 
less, it was one of the motivations for the presented study of analyticity under the branch 
cut in Fig. 2. The “extended analyticity” proposed in Sect. 4 is the statement that the 
nearest singularity under this branch cut is identical to the Yang-Lee singularity in Fig. 1. 
In this terms, the “true spinodal” is the Yang-Lee singularity (viewed from the shadow 
side, as in Fig. 4). 

Our analysis concerned explicitely with the D = 2 system. In this case the Yang-Lee 
singularity turns out to be located relatively far under the branch cut. The phase of ^sp 
in Fig. 2 (which is related to the angular size of the “shadow domain” in Fig. 4) is Stt/S. 
However, on a formal level, it is not difficult to extend the analysis of Sect. 3 to any 
1 < T> < 4. In this generic case the angular size 9 of analogous shadow domain can be 
expressed in terms of the standard Ising model critical exponents as 



1 

2 — a — (3 


( 9 . 1 ) 


It follows from what we know about the T>-dependence of the Ising critical exponents 
(see e.g. Ref. [^) that when D increases from 2 to 4 the angular size of the “shadow 
domain” decreases, and goes to zero as the critical exponents approach their classical 
values at D = 4. A simple estimate shows that at D = 3 this angular size is ~ 0.047r. 
If similar “extended analyticity” assumption could be validated in this case, one should 
expect the effect of the Yang-Lee singularity being quite prominent at real negative H, 
at least in the critical domain T ^ Tc- Unfortunately, at the moment we do not see any 
practical way to check this assumption away from the D = 2 case. 

Even in the D = 2 case many issues call for rehnement. First, although the extended 
dispersion relation did check in Sect. 9 rather convincingly, with present data about the dis¬ 
continuity A the agreement in Table 4 is not too precise, and hence very weak singularities 
in the shadow domain are not excluded. Much better knowledge about the “hne-structure” 
of the Yang-Lee singularity (in particular, better estimates of the coefficients of the ex¬ 
pansions ( |3.41fc|) , ( |3.41cj) ) would be usefull in order to establish the extended analyticity 
“beyond any reasonable doubt”. We hope to make progress in this direction in the fu¬ 
ture. We also want to stress here that even that would not settle the question of analytic 
properties of the scaling function completely. What one should expect when going under 
the Yang-Lee branch cut in Fig. 4 from the shadow domain? This and related questions 
require at least some understanding of the physics governing the “shadowy side” of the 
Yang-Lee edge singularity. 
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Let us also mention some possible further developments, not related directly to the 
free energy, which we hope to address in the future: 

(t) Analytic properties of the correlation length (dehned as inverse mass Mi), can 
be studied along the same lines. We expect associated scaling function to enjoy similar 
extended analyticity. 

(m) Detailed quantitative study of the IFT (|1.1|) , viewed as the particle theory, is 
possible using TFFSA. This is a beautifull model of quark conhnement, simple but rich in 
phenomena . 


{iii) Physical interpretation of the IFT with pure imaginary h is a particularly in¬ 
triguing question. Even the qualitative pattern of the higher hnite-size energy levels Ei[R) 
in this case is yet to be understood. 


Appendix A. 


Here we present a simple derivation of the hnite-size matrix elements (|2.12|) using the 
symmetries of the “doubled” Ising held theory with h = 0. This derivation is one of the 
results of |^, where these symmetries are exploited to a greater extent (in particular, the 
nonlinear diherential equations of 0 are shown to be a consequence of these symmetries). 

Consider two copies of the free-fermion system (|2.2p . We will use the subscripts a 
and b to distinguish between the helds corresponding to these copies. Thus, (V^a,'0a) and 
'ijjb) will stand for the two species of Majorana fermi helds, aa and Uf, will denote 
the associated spin helds, etc. We will also use the notations a^, a and 6^, b for the 
fermion creation and anihilation operators associated with the two copies; these operators 
are assumed to satisfy the canonical anticommutators (|2.6| ) between themselves, and a^, a 
anticommute with &!, b. Of course, this doubled Majorana fermion system is equivalent to 
a single copy of a free Dirac fermion. 

We will study this system in a hnite-size geometry, with the spatial coordinate x 
compactihed on a circle of circumference R. Let us call Tidiag the space of states generated 
by the creation operators from the “diagonal” vacua | 0)ns,ns =| 0)NS,a® I 0)NS,f) 

and I 0 )r,r =1 0)R,a{8) | 0)R,f,. Consider the operator 


Y = i- 

27r 


-i? r 


XfYl — 

ll^ad'lpb + '>Pai>b 


dx 


(A.l) 


acting in the space 'Hdiag- It is easy to check that it commutes with the Hamiltonian of 
this doubled system, and anihilates both the above “diagonal” vacua | 0)ns,ns and | 0)r^r. 
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This operator is just one of the generators of the inhnite-dimensional algebra SL{2) of the 
local integrals of motion of the free Dirac held theory [^; in fact, for our present limited 
purpose no knowledge about the other generators is needed. It is straightforward to derive 
the commutators 


[Y, iia{x)nh{x)] = aa{x) dah{x) - daa{x) ah{x ), 

[Y, aa{x)iJib{x)] = -i {iia{x) dab(x) - dna{x) (Jb{x)) , 


where again d stands for — idy), and n{x) are the dual spin held 


(A.2a) 

(A.26) 

It is also 


straightforward to compute the action of this operator on arbitrary states containing a 
and b particles. For instance, for the excitations over the vacuum | one obtains 

m 


Y = [alb^ + b^a^] , 


(A.3) 


nez 


k G Z + i) holds for the 


where dn are dehned in ( p.l3|) . Similar expression (with n 
action of Y in the NS-NS sector. We are going to show that the general structure of 
the matrix elements ( p.l2|) (with the exception of the explicit expressions ( p.l9|) for the 
“leg factors” g and g, which will be derived in the second part of this section) is a simple 
consequence of the relations (A.2) and ( A-3|) . 

To save on writing, we will restrict attention to the case m > 0, and explicitely derive 


only the particular case A = 0 of the expression ( 2.12 ); the analysis can be extended to a 


f NS(0 

ct(0,0) 

ni,. 

• nN)R 

for even 

N] 

(A.4) 





i NS(0 

^(0,0) 

ni, . 

nN)R 

for odd 

Nj 

; (0 1 Y, Ha 

{x)gb{x] 

■ 

^riN 

0)r,r . 


(A.5) 


other cases in a straightforward way. We will use the notations 
Mo,w(ni,..., un) 

If N is even, take the identity 


which follows from ns,ns(0 | Y = 0. If Y is odd, similar identity, with ga{x)fib{x) replaced 
by cra{x)fib{x), should be used. The matrix element in the r.h.s. of (|A.5|) can be written 
as a sum of the matrix element of the commutator (|A-2aD and the matrix element of 
fJ'aix)gb{x)Y . Using (|A.3|) to compute the latter, one obtains the relation 

N 

Afo,oMo,Ar(ni,---,nAr) == 

I ^ (-1)^“^ Mo,i(nj) Mo,Ar-i(ni, • • •, rij-i, rij+i, • • •, Un) , 
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where / = i if is even, and / = 1 if is odd. This relation allows one to obtain 
recurrently all the matrix elements (|A.4|) in terms of Mo,o and Mo,i(n). It is not difficult 
to show that the solution of this recurrent relation has the form 


N 


^0,Ar(ni, • • • , nw) = Mo,o I 5'(^ni) ) fN{.0n^, ■ ■ ■ ,6nj^) , (A.7) 




where [A^/2] denotes the integer part of N/2, 


g{6n) = Mo,i(n)/Mo,o, 


(A.8) 


and 


N 


/jv(6'i,---,6'jv) = JJ tanh 


i<j 


Om - Ouj 


(A.9) 


This follows from the well known identities ||^ that the products (|A.9|) obey, 

N 


fo,N+i{ 0 , 0 i,-■ ■ , 0 n) ='^{-y ^ tanh 

i=i 


9 - 9 , 


fo,N—ii9ij ■ ■ ■ 5 ^j—it 9j+if ■ • • , 9j\i') 


(A.IO) 

If one expands ( |AT0|) around 9 = +cxd, the term ~ e ^ yields exactly (|A.6|) . 

It remains to determine the factors Mq^q and g{ 9 ) in ( |A.7| ). The hrst is known from 
Ref. 1 ^, see Eqs. ( |2.1ti|) -( p.l7| ). To hnd the factors g{9n) =ns ( 0 |;a( 0 )|n)R/Ns( 0 |cT( 0 )| 0 )R, 
let us return to the single Majorana fermion system (|2.2|), and consider the matrix element 


T(x,y) =NS (0 I /r(0,0) V^(x,y) | 0)r/ns(0 | cr(0, 0) | 0)r, 


(A.ll) 


and similar matrix element T(x, y) defined as in Eq. ( A.ll ) with '0(x, y) replaced by '0(x, y). 
Here x, y are the cartesian coordinates on the euclidean space-time cilinder, x being the 
periodic direction. The functions T and T have the following properties: 
i) They satisfy the Majorana held equations 


a9<(x,y) = !^f(x,y), 54'(x,y) = -^^^(x.y); 


(A. 12) 


ii) T(x, y) is an analytic function for all real x, y except at the points (x, y) = (0(mod R), 0), 
and it obeys the (anti)periodicity conditions 


4'(x + R,y) = 


T(x, y) for y<0 
-\l/(x,y) for y>0 


(A. 13) 
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The function T obeys the same conditions. 

in) Their singularities at (x, y) = (0, 0) are such that 


lim a/^T(x, y) = e , lim \/^T(x, y) = e 4 " , (A.14) 

x,y^0 x,y^0 

where ^ = x + fy, z = y. — iy. Also, the functions T(x, y) and T(x, y) are bounded when 
|y| ^ cx). 

It is possible to show that the above conditions i)-iii) dehne the functions uniquely. 
The solution to these conditions has the form 


^'(x,y) 

^'(x,y) 


n6Z 

E 

fc62+| 


TT 


R cosh 9r 


TT 

R cosh 6*1 


e 2 " g{6n) e 


my cosh On—imx sinh^^i 


e 2 g{ek) e 


—my cosh^fc+zmx sinh^fc 


for y < 0 , (A.15a) 
for y > 0 , (A.156) 


where 6 *„ and 9k are the solutions of the Eqs. (|2.13| ), and the functions g{9) and g{9) 
are dehned through the integrals (|2.18|) . The function T(x, y) is given by the same sums 
( |A.15a|) , (|A.156| ), but with the factor ei replaced by e“5. It is clear that the sums (|A.15a|) , 
( |A.156|) satisfy ([A.12|) -( |A.14D , and one only needs to check that the analytic continuation 
of (A-15a|) to y > 0 coincides with ( |A.156|) , and vice versa. This is done using standard 


tricks with contour deformations. One can write (A-15a|) , (|A.156| ) in terms of the contour 
integrals 


4.(x,y)= ^ fory<0,(A.16<.) 

<P(x,y) = *y?S^ fory>0, (A.16!,) 

where G{9) = exp(ft;( 6 *)), with k,{9) dehned as the integral (|2.19| ), G{9)~^ stands for 1/G(6*) 
(not for the inverse function), and the integration is performed along the contours C_, which 
runs from —00 to +cxo just below the real axis, and C+, which returns back to —00 just 
above the real axis. The function G{9) is meromorphic, and both G{9) and G{9)~^ are 
analytic in the strip —tt < Q'm 9 < n. It is easy to verify that 

pTnR cosh 9 _ i 

0(0 + ^^/2)G{9 - Z./ 2 ) = . (A.17) 
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If, for instance, 0 < x < R, then the contour C_ in ( A.16a ) and ( A.166 ) can be shifted 
downward, 9 ^ 9 — irr j2, while the contour C-t- admits upward shift 9 9 R it: j2. After 

this deformation the integral (|A.16a|) assumes the form 


^'(x, y) = \pnm 


d9 

-r^ 2TI 


G{9 in j2') e 2 4 slnh d—mx cosh 9 


_ g —mi? cosh 0 

G{9 + m/2) 


^imy sinh 6-\-mx cosh 6 


1 


^mi? cosh 9 


(A. 18) 


As long as x remains within the above interval, this integral now converges at all y, both 
negative and positive; it thus dehnes the analytic continuation of (|A.16a|) to all y. The 
same contour deformations applied to ( A.166 ) yields the integral which coincides with 
( |A.18|) in virtue of ( |A.17| ). 

In fact, it is not difficult to hnd in a similar way a complete basis of functions \l/(x, y) 
which satisfy (|A.12|) - (|A.14|) , but are allowed to grow at |y| —cx); this basis then can be 
used to derive the full set of matrix elements (|2.12|) directly. We will discuss this topic in 
greater detail elsewhere. 


Appendix B. 

In this section we describe some properties of the energy spectrum of the Hamilto¬ 
nian ( p.3|) in the low-T domain m > 0. In subsection BI below we warm up with the 
derivation of the small h expansion ( |5.2|) of the meson masses, and in BII we discuss 
the relation of the near-intersection pattern in Fig. 6a to the characteristics of the “false 
vacuum” resonance. 


BI. As is well known, in the low-T regime the interaction term h f a(x) d?x in (2.1) 
gives rise to a conhning interaction between the fermions of (|2.2| ) (the “quarks”). As the 
result, the particle spectrum of the IFT ( prT|) with h Q consists of bound states of the 
quarks - the “mesons”. If /i is sufficiently small, the lower part of the meson spectrum can 
be studied within the two-quark approximation. One looks for the meson state (in its rest 
frame) in the form 

r ^ 

Stt 


T) = 


T(p) ajalp 


0 ), 


(B.l) 


and ignores all multiquark components. Here aj, are the quark creation operators 
normalized according to the canonical anticommutators {a^, a^,} = 27r 5{p — p'), and 
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^(p) is the momentum-space wave function, which is assumed to be antisymmetric, i.e. 

Within this approximation the eigenvalue problem for the Hamilto¬ 
nian (|2.3|) (with R = oo) leads to the Bethe-Salpeter equation 


{2u{p) - E) ^{p) = 


3 roo 


m 


2 7-00 


u{p) + u{p')\^ ^ 1 pp' 


p — p 


2 uj{p) w(p') 


(B.2) 

where the principal value of the integral in the right-hand side is understood, E stands for 
the meson’s rest energy, and the notations ijo{p) = and 


, 2ah \ 


(B.3) 


are used. The easiest way to derive this equation is to use the hnite-size matrix ele¬ 
ments ( ^.121 ) with N = M = 2 and take the limit R ^ oo. 

The two-particle approximation is justihed if h is small and E is sufficiently close 
to 2m. In this case one can consistently treat the momenta p, p' in (p.2|) as small as 
compaired to m. Making a rescaling 


p = {mt) k , p' = {mt) k ', 


(B.4) 


and expanding the operators in (B2) in the powers of t, one has 

T((mt) k) 


f2 fA 

-e + k^ - —k'^ R — k^ ^ - 

4 8 




? / / . , dk' 

T((mt)/c ) — , 


where e is dehned as 


E — 2m = {mt^) e . 

It is useful to write down the conhguration-space form of the Eq. (p.5|) . 


(f (f 

^ ~ ~ T ~ s’^ 


where the conhguration-space wave function 

r dk 


i,(X) = 




(B.5) 


(B.6) 


V>(V), (B.7) 


i,(X) 


AkX 


2ti 


T((mt) k) 


(B.8) 
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is written in terms of the rescaled coordinates, 


X = (mt)(xi - X2), 


(B.9) 


xi,X 2 being the positions of the qnarks. This function also is antisymmetric, 


i,(x) = -V>(-X) 


(B.IO) 


In Eqs. ([b.5| ) and (|B.7|) the omitted terms ■ ■ ■ are ~ or smaller. 

In the leading order in the Eq. ( p.7|) is just the Schroedinger equation with the 
linear potential |X|. It leads to the assymptotic formula (|5.1|) for the meson masses (hrst 
obtained in |^]). One can calculate further corrections by solving ( B.7 ) perturbatively 
in Let A{X) be a solution of the Airy equation 


A - 




dX^ 


A{X) = 0. 


(B.ll) 


It is straightforward to check that the function 


f(x)=a(x)-L 


AX A{X) + X^ A'{X) 


A 


7 200 / ^ ^ V 35 35 ' ^ ^ 


(B.12) 




(where prime denotes the derivative) is the perturbative solution of the equation 

d^ e d^ d^ 


X 


and that the function 


dX^ 4 dX^ 8 dX^ 


+ 0{A 


F{X) = 0, 


i;{X)=sign{X)Fi\X\-e) 

solves (in the same perturbative sense) the equation (IB.711, provided 


(B.13) 


(B.14) 


F{-e) = 0(e ); 

F(-e) + ^F"(-e)=0(t4); 








F{-e) + ^ F"(-6) + ^ ^ F'{-e) = 0{t^). 


(B.15a) 

(B.156) 

(B.15c) 
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Now, let Ai(X) and Bi(A) be the two canonical solntions of the Airy eqnation ( [B.ll ), and 
Fa{X) and Fb{X) - the associated pertnrbative solntions ( |B.12|) of ([B.13| ). A straightfor¬ 
ward calcnlation shows that the linear combination 


with 


20 



F(A) = 

Fa{X) + 

a(e) = 

Ai(-6) 

r 

^ 20 

Bi(-e) 

Bi'(-6) 

Ai'(-6)\ 

(2e^ 57 

Bi(-e) 

Ai(-6) ; 

iy 35 14 


Bb(-6) Ab(-6) 

Bi(-e) Ai(-e) 


1 Bi' 


20 Bi(- 


+ 0{t^ 


(B.16) 


(B.17) 


satishes the conditions (|B.15a|) -( [B.15c|) , i.e. ([B.14|) with this F{X) solves ([B.7|) . 


In the bonnd-state problem one is interested in the normalizable solntions of (|B.7|) , 

iZc \^{X)\'^dX < oo, i.e. 


ip{X) —>0 as 


|A| 


CX) . 


(B.18) 


Threfore, the bonnd-state eigenvalnes of the energy parameter e in (|B.7|) coincide with the 
zeroes of the amplitnde cr(e) in (|B.16|). From (|B.17|) one hnds for these zeroes, 


ei = 


20 


280 
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y ^? ) + 0(i'=' 


(B.19) 


where —Zi, i = 1,2,3,, are the zeroes of the Airy function, Ki{—Zi) = 0. 

It is straightforward to continue this perturbative expansion to higher orders in 
However, without proper modihcations this does not make much sense because, starting 
with the order in (|B.7|) , this expansion exceeds the accuracy of the two-particle approxi¬ 


mation. The exact meson state | T) contains four-quark, six-quark, and higher components 
which are neglected in ( [B.l|) . One can check that these multi-quark components lead to 
certain corrections to the r.h.s. of (|B.2|) which start with the terms ~ Therefore, 
strictly speaking, even the terms ~ t'^ in (|B.7|) and in ( p.l9|) go beyond the accuracy of the 
two-quark approximation. However, it is possible to argue that the whole of the leading 
(~ t®) correction to (p.2|) due to the multi-quark contributions has the form of the quark 
mass renormalization ( |5.8| ), i.e. it ammounts to adding a momentum-independent term 
Q 2 I 2 to the l.h.s. of ([B.7|) . In turn, this leads to the shift -|- q2/2 in (|B.19|) . 

Using 


one arrives at the expansion (b.2|). 
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BII. In the low-T domain m > 0 the gronnd state of the IFT with R = oo and h = 0 is 
two-fold degenerate, with the two vacnnm states | 0)± corresponding to the oposite valnes 
of the spontaneons magnetization, 


± 


(0 I cr{x) I 0)± = ± cr. 


(B.20) 


Adding the interaction term h J <j{x) dPx with small positive h lowers the energy density 
associated with the state | 0)_ and rises the energy density of | 0)_(_. That is, the state | 
0)_ becomes the trne vacnnm of the perturbed system, while the state | 0)-|- looses its 
stability and becomes the “false vacuum” - the global resonance state with the complex 
energy (|3.16|) . 

If h is small, the “false vacuum” decay rate can be computed using the saddle-point 
analysis in the Euclidean space-time [p9|,^,|iI|,pT| . This analysis suggests the following 


qualitative picture of the decay process in real time. The quantum mechanical tunneling 
process results in the simultaneous formation of two quarks separated by the distance Rd = 
2m/AF ^ m/ah, which then speed away from each other, the expanding region between 
the quarks being the nucleus of the true vacuum. The distance Rd is the size of the “critical 
droplet” in the droplet model calculation [^,|^. As long as the separation Ri ^2 between 
the quarks remains close to Rd, the state can be described in term of non-relativistic 
quantum mechanics of the two-quark system with the repulsive linear potential —2ha 
If |i?i ,2 — Rd\ ^ Rdi a generic state of such two-quark system is described by the wave 
function 

^{Y) ~ A+(y - e') + Sres{E) A_(Y - e'), (B.21) 

where Y = {mt) {Rd — Ri, 2 ) (the parameter t is dehned in (|B.3| )), {mt^) e' = E — Ej-es 
is suitably normalized energy of the system over the energy Eres of the “false vacuum”, 
and A± are the solutions of the Airy equation ( [B.llj) , 


A±(y) = Ai(y) ±fBi(y). 


(B.22) 


These are the “running wave” solutions, the wave A-{Y) being the state of the quarks 
speeding away from each other. The “scattering amplitude” Sres in ([B.21| ) should be de¬ 
termined (in principle) by solving the problem in the domain where the quarks are deep 
inside the classically unaccesible region Ri ^2 < Rd- In this region one expects the rela¬ 
tivistic effects, in particular the multi-quark components of the state, to play a signihcant 
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role. The solution taking into account these effects is not yet availablel, and the exact 
form of Sres is not known. However, even in the absence of the explicit solution, one can 
predict that the amplitude Sres{E), as the function of the total energy E of this state, must 
exhibit a resonance pole at the complex energy (|3.16 ). If E is close to the resonance energy 
(given by the real part of (|3.16|) ), this amplitude can be approximated by the Breit-Wigner 
formula 

E/R- AE + iT 


Sres{E) 


(B.23) 


E/R- /\E -iV ’ 

where R is the spatial size of the system, and AT = E^es/R — 2 d h + 0{h^) is the same 
as in (O). With this approximation for the “resonance scattering amplitude” Sres{E) 
in (p.21|) , and with the relativistic corrections described in BI taken into account, the 
wave function (p.21|) can be written as 


■0(1") = const. - e') + ares(e) Eb{Y - e')^ , 


where 


(0 = 


r(m, h) R 
mt'^ e’ 


(B.24) 


(B.25) 


and Ea and Eb are the same perturbative solutions of (|B.13|) as in ([B.16|) . 

Now let us come back to the problem of the meson states, this time assuming a finite- 
size geometry of the system, with the spatial coordinate x compactified on a circle of finite 
circumference R. If h is sufficiently small, almost all the analysis of BI remains valid. 


with the exception of the bound-state condition ([B.18 ), which has to be slightly modified. 
When the meson (i.e. the two-quark system with the confining interaction generated by 
the magnetic field h) is put on a finite-size circle, its energy levels should be determined 
from some sort of periodicity condition imposed on the state wave function, rather then by 
( |B.18|) . Again, unfortunately the two-quark approximation in general breaks down in the 
classically unaccessible domain |A| e (here |A| is the rescaled quarks separation (|B.9| ), 
and e relates to the energy E as in ([B.6|) ), where multi-quark components of the wave 
function become significant. Although the solution for the full wave function in this domain 
is not known, one can note that this problem is closely related to the above problem of 


^ An attempt to apply the two-particle approximation in this region was made in an interestiong 
recent paper We comment upon the result of this work at the end of Sect. 7.2. 
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determining the “scattering amplitude” Sres in ( B.21 ), and the form ( B.25 ) of the wave 
function still applies if one makes identifications 


e' = e — mt {R — Rq 


Y -e' = X 


(B.26) 


which follow from the periodicity of the problem. With this, by demanding agreement 
between (|B.16D and ([B.24|) , one arrives at the equation 


a{e) = cTres(e -mt{R- Rd )), 


(B.27) 


which determines the R dependence of the hnite-size energy levels Ei{R) — Eq{R) = 2m + 
mt^ ei{R) associated with the meson states; here ei{R) are the solutions of the equation 

( 133 ) ■ 

The equation ( p.27|) is derived within the two-quark approximation which is justified 
if h is small. It applies when R is sufficiently greater then Rd (otherwise the domain of 
validity of ( |B.24| ) does not exist). It is not difficult to see that in this domain ([B.27|) leads 
to a pattern of the hnite-size energies very similar to what is shown in the Fig. 6a. Indeed, 
for small h the resonance width F is expected to be exponentially small (see Eq. ( 3.18D ), 
and hence the function ares in the r.h.s. of (|B.27| ) is very small everywhere exept when e is 
close to its resonance pole eres{R) = nrt (R — Rd). Therefore, Eq. (|B.27|) can hold either by 
virtue of e being close to one of the zeroes of the function a(e), which leads to E{R) — Eq{R) 
close to one of the meson masses M^, or by virtue of e being close to eres(-R), leading to 
E{R) — Eo{R) close to AE R. When the both conditions are met, which happens when R 
approaches one of the “near-intersection” points Ri,i+i = Mi/AE, one can use the linear 
approximation a(e) ~ (e — e^) a'{ei), which results in a quadratic equation for e{R). This 
way one arrives at Eq. (|5.13|) for the separations Ei+i{R) — Ei{R). 
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